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Preface 


Electromagnetism is fundamental to the whole of electrical and electronic engineering. It provides the 
basis for understanding the uses of electricity and for the design of the whole spectrum of devices 
from the largest turbo-alternators to the smallest microcircuits. This subject is a vital part of the 
education of electronic engineers. Without it they are limited to understanding electronic circuits in 
terms of the idealizations of circuit theory. 


The book is, first and foremost, about electromagnetism, and any book which covers this subject must 
deal with its various laws. But you can choose different ways of entering its description and still, in 
the end, cover the same ground. I have chosen a conventional sequence of presentation, beginning 
with electrostatics, then moving to current electricity, the magnetic effects of currents, 
electromagnetic induction and electromagnetic waves. This seems to me to be the most logical 
approach. 


Authors differ in the significance they ascribe to the four field vectors E, D, B and H. I find it 
simplest to regard E and B as ‘physical’ quantities because they are directly related to forces on 
electric charges, and D and H as useful inventions which make it easier to solve problems involving 
material media. For this reason the introduction of D and H is deferred until the points at which they 
are needed for this purpose. 


Secondly, this is a book for those whose main interest is in electronics. The restricted space available 
meant that decisions had to be taken about what to include or omit. Where topics, such as the force on 
a charged particle moving in vacuum or an iron surface in a magnetic field, have been omitted, it is 
because they are of marginal importance for most electronic engineers. I have also omitted the chapter 
on radio-frequency interference which appeared in the second edition despite its practical importance. 


Thirdly, I have written a book for engineers. On the whole engineers take the laws of physics as 
given. Their task is to apply them to the practical problems they meet in their work. For this reason I 
have chosen to introduce the laws with demonstrations of plausibility rather than formal proofs. It 
seems to me that engineers understand things best from practical examples rather than abstract 
mathematics. I have found from experience that few textbooks on electromagnetism are much help 
when it comes to applying the subject, so here I have tried to make good that deficiency both by 
emphasizing the strategies of problem-solving and the range of techniques available. A companion 
volume is planned to provide worked examples. 


Most university engineering students already have some familiarity with the fundamentals of 
electricity and magnetism from their school physics courses. This book is designed to build on that 
foundation by providing a systematic treatment of a subject which may previously have been 
encountered as a set of experimental phenomena with no clear links between them. Those who have 
not studied the subject before, or who feel a need to revise the basic ideas, should consult the 
elementary texts listed in the Bibliography. 
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The mathematical techniques used in this book are all covered either at A-level or during the first year 
at university. They include calculus, coordinate geometry and vector algebra, including the use of dot 
and cross products. Vector notation makes it possible to state the laws of electromagnetism in concise 
general forms. This advantage seems to me to outweigh the possible disadvantage of its relative 
unfamiliarity. I have introduced the notation of vector calculus in order to provide students with a 
basis for understanding more advanced texts which deal with electromagnetic waves. No attempt is 
made here to apply the methods of vector calculus because the emphasis is on practical problem- 
solving and acquiring insight and not on the application of advanced mathematics. 


I am indebted for my understanding of this subject to many people, teachers, authors and colleagues, 
but I feel a particular debt to my father who taught me the value of thinking about problems ‘from 
first principles’. His own book, The Electromagnetic Field in its Engineering Aspects (2nd edn, 
Longman, 1967) is a much more profound treatment than I have been able to attempt, and is well 
worth consulting. 


I should like to record my gratitude to my editors, Professors Bloodworth and Dorey, of the white and 
red roses, to Tony Compton and my colleague David Bradley, all of whom read the draft of the first 
edition and offered many helpful suggestions. I am also indebted to Professor Freeman of Imperial 
College and Professor Sykulski of the University of Southampton for pointing out mistakes in my 
discussion of energy methods in the first edition. 


Finally, I now realize why authors acknowledge the support and forbearance of their wives and 
families through the months of burning the midnight oil, and I am most happy to acknowledge my 
debt there also. 


Richard Carter 
Lancaster 2009 
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1. Electrostatics in free space 


Objectives 


e To show how the idea of the electric field is based on the inverse square law of force between 
two electric charges. 

e To explain the principle of superposition and the circumstances in which it can be applied. 

e To explain the concept of the flux of an electric field. 

e To introduce Gauss’ theorem and to show how it can be applied to those cases where the 
symmetry of the problem makes it possible. 

e To derive the differential form of Gauss’ theorem. 

e To introduce the concept of electrostatic potential difference and to show how to calculate it 
from a given electric field distribution. 

e To explain the idea of the gradient of the potential and to show how it can be used to calculate 
the electric field from a given potential distribution. 

e To show how simple problems involving electrodes with applied potentials can be solved 
using Gauss’ theorem, the principle of superposition and the method of images. 

e To introduce the Laplace and Poisson equations. 

e To show how the finite difference method can be used to find the solution to Laplace’s 
equation for simple two-dimensional problems. 
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1.1 The inverse square law of force between two electric charges 


The idea that electric charges exert forces on each other needs no introduction to anyone who has ever 
drawn a comb through his or her hair and used it to pick up small pieces of paper. The existence of 
electric charges and of the forces between them underlies every kind of electrical or electronic device. 
For the present we shall concentrate on the forces between charges which are at rest and on the force 
exerted on a moving charge by other charges which are at rest. The question of the forces between 
moving charges, which is a little more difficult, is dealt with in Chapter 4. 


The science of phenomena involving stationary electric charges, known as electrostatics, finds many 
applications in electronics, including the calculation of capacitance and the theory of every type of 
active electronic device. Electrostatic phenomena are put to work in electrostatic copiers and paint 
sprays. They can also be a considerable nuisance, leading to explosions in oil tankers and the need for 
special precautions when handling metal-oxide semiconductor integrated circuits. 


The starting point for the discussion of electrostatics is the experimentally determined law of force 
between two concentrated charges. This law, first established by Coulomb (1785), is that the force is 
proportional to the product of the magnitudes of the charges and inversely proportional to the square 
of the distance between them. In the shorthand of mathematics the law may be written 


F= LQ: 5 (1.1) 
Are yr 

where Q, and Q) are the magnitudes of the two charges and r is 

the distance between them, as shown in the figure on the right. 

Now force is a vector quantity, so Equation (1.1) includes the 


unit vector F which is directed from Q; towards Q, and the 


equation gives the force exerted on Q> by Q;. The force exerted 
on Q, by Q> is equal and opposite, as required by Newton’s 
third law of motion. 


Examination of Equation (1.1) shows that it includes the effect of the polarity of the charges correctly, 
so that like charges repel each other while unlike charges attract. The symbol é denotes the primary 
electric constant; its value depends upon the system of units being used. In this book SI units are 
used throughout, as is now the almost universal practice of engineers. In this system of units é is 
measured in Farads per metre, and its experimental value is 8,854 x 10!” F m’'; the SI unit of charge is 
the coulomb (C). 
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Electric charge on a macroscopic scale is the result of the accumulation of large numbers of atomic 
charges each having magnitude 1.602 x 10°'° C. These charges may be positive or negative, protons 
being positively charged and electrons negatively. In nearly all problems in electronics the electrons 
are movable charges while the protons remain fixed in the crystal lattices of solid conductors or 
insulators. The exceptions to this occur in conduction in liquids and gases, where positive ions may 
contribute to the electric current. 


1.2 The electric field 


Although Equation (1.1) is fundamental to the theory of electrostatics it is seldom, if ever, used 
directly. The reason for this is that we are usually interested in effects involving large numbers of 
charges, so that the use of Equation (1.1) would require some sort of summation over the (vector) 
forces on a charge produced by every other charge. This is not normally easy to do and, as we shall 
see later, the distribution of charges is not always known, though it can be calculated if necessary. 
Equation (1.1) can be divided into two parts by the introduction of a new vector E, so that 


p-—2 3 


= 1.2 
Aner’ (1.2) 


and F=Q.E (1.3) 


The vector E is known as the electric field, and is measured in volts per metre in SI units. The step of 
introducing E is important because it separates the source of the electric force (Q;) from its effect on 
the charge Q2. The question of whether the electric field has a real existence or not is one which we 
can leave to the philosophers of science; its importance to engineers is that it is an effective tool for 
solving problems. 


The electric field is often represented by diagrams like Fig. 1.1 in which the lines, referred to as ‘lines 
of force’, show the direction of E. The arrowheads show the direction of the force which would act on 
a positive charge placed in the field. The spacing of the lines of force is close where the field is strong 
and wider where it is weak. This kind of diagram is a useful aid to thought about electric fields, so it is 
well worth while becoming proficient in sketching the field patterns associated with different 
arrangements of charges. We shall return to this point later, when discussing electric fields in the 
presence of conducting materials. 
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Fig.1.1 The electric field of a point charge can be represented diagrammatically by lines of force. The 
figure should really be three-dimensional, with the lines distributed evenly in all directions. 


In order to move from the idea of the force acting between two point charges to that acting on a charge due to 
a whole assembly of other charges it is necessary to invoke the principle of superposition. This principle 
applies to any linear system, that is, one in which the response of the system is directly proportional to the 
stimulus producing it. The principle states that the response of the system to a set of stimuli applied 
simultaneously is equal to the sum of the responses produced when the stimuli are applied separately. 


Equation (1.2) shows that the electric field in the absence of material media (‘in free space’) is proportional 
to the charge producing it, so the field produced by an assembly of charges is the vector sum of the fields 
due to the individual charges. The principle of superposition is very valuable because it allows us to tackle 
complicated problems by treating them as the sums of simpler problems. It is important to remember that 
the principle can be applied only to linear systems. The response of some materials to electric fields is non- 
linear and the use of the principle is not valid in problems involving them. 
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Before discussing ways of calculating the electric field it is worth noting why we might wish to do it. 
The information might be needed to calculate: 

e the forces on charges; 

e the conditions under which voltage breakdown might occur; 

e capacitance; 

e the electrostatic forces on material media. 
The last of these is put to use in electrostatic loudspeakers, copiers, ink-jet printers and paint sprays 
and is of growing importance in micro-mechanical devices. 


1.3 Gauss' theorem 


Figure 1.1 shows electric field lines radiating from a charge in much the same way that flow lines in 
an incompressible fluid radiate from a source (such as the end of a thin pipe) immersed in a large 
volume of fluid (shown in Fig. 1.2). Now in the fluid the volume flow rate across a control surface 
such as S, which encloses the end of the pipe, must be independent of the surface chosen and equal to 
the flow rate down the pipe, that is, to the strength of the source. 


Fig. 1.2 The flow of an incompressible fluid from the end of a thin pipe is analogous to the electric 
field of a point charge. It is necessary for the end of the pipe to be well away from the surface of the 
fluid and the walls of the containing vessel 


To apply this idea to the electric field it is necessary to define the 

equivalent of the flow rate which is known as the electric flux. The figure dA 
on the right shows a small element of surface of area d4 and the local 
direction of the electric field E. The flux of E through dA is defined as the E 
product of the area with the normal component of E. This can be written 
very neatly using vector notation by defining a vector dA normal to the 
surface element. The flux of E through dA is then just E-dA = E dd cos 8. 
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Now consider the total flux coming from a point charge. The simplest choice of control surface 
(usually called a Gaussian surface in this context) is a sphere concentric with the charge. Equation 
(1.2) shows that E is always normal to the surface of the sphere and its magnitude is constant there. 
This makes the calculation of the flux of E out of the sphere easy - it is just the product of the 
magnitude of E with the surface area of the sphere: 


flux of E= Fi Arr’ == (1.4) 


Thus the flux of E out of the sphere is independent of the radius of the sphere and depends only on the 
charge enclosed within it. It can be shown that this result is true for any shape of surface and, by using 
the principle of superposition, for any grouping of charges enclosed. The result may be stated in 
words: 


The flux of E out of any closed surface in free space is equal to the charge enclosed by the surface 
divided by €o. 


This is known as Gauss’ theorem. It can also be written, using the notation of mathematics, as 


ff Eda=—Iff par (1.5) 
S OV 


where S is a closed surface enclosing the volume V and p is the charge density within it. Equation 
(1.5) looks fearsome but, in fact, it is possible to apply it directly only in three cases whose symmetry 
allows the integrals to be evaluated. Those cases are: 

e Parallel planes 

e Concentric cylinders 

e Concentric spheres 


To show how this is done let us consider the case of a long, straight, rod of radius a carrying a 
uniform charge q per unit length. From the symmetry of the problem we can assume that E is 
everywhere directed radially outwards and that the magnitude of E depends only on the distance from 
the axis. This is not valid near the ends of the rod but the problem can be solved in this way only if 
this assumption is made. The next step is to define the Gaussian surface to be used. This is chosen to 
be a cylinder of radius r and unit length concentric with the charged cylinder with ends which are flat 
and perpendicular to the axis as shown in Fig. 1.3 
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Fig. 1.3 A Gaussian surface for calculating the electric field strength around a charged rod. 


On the curved part of the Gaussian surface E has constant magnitude and is everywhere perpendicular 
to the surface. The flux of E out of this part of the surface is therefore equal to the product of FE and 
the area of the curved surface. On the ends of the cylinder E is not constant but, since it is always 
parallel to the surface, the flux of E out of the ends of the cylinder is zero. Finally, since the Gaussian 
surface is of unit length it encloses charge g. Therefore, from Gauss’ theorem, 


2arE =q/é&, 
or E=q/2ne yr (1.6) 
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1.4 The differential form of Gauss’ theorem 


Only a limited range of problems can be solved by the direct use of Equation (1.5). Another form, 
which is obtained by applying it to a small volume element, enables us to solve a much wider range of 
problems. Figure 1.4 shows such a volume element in Cartesian coordinates. 


x 


Fig. 1.4 The elementary Gaussian surface used to derive the differential form of Gauss’ theorem. 


To calculate the net flux out of the element, consider first the two shaded faces A and B which are 
perpendicular to the x-axis. The only component of E which contributes to the flux through these 
faces is F.. If the component of E is £, on A then, in general, the x component of E on B can be 


written EL + (OE el ax) dx . Provided that the dimensions of the element are small enough we can 


assume that these components are constant on the surfaces A and B. The flux of E out of the volume 
through the faces A and B is then 


oH iE 
ae dx i= aye 
) Ox 


IX 


—E, dydz+ G + 


The same argument can be used for the other two directions in space, with the result that the net flux 
of E out of the element is 


OE 
d® , = TE, Oi Os dx dy dz (1.7) 
Ox Ov 


Now if p is the local charge density, which may be assumed to be constant if the volume element is 
small enough, the charge enclosed in the volume is 


dO = pdxdydz (1.8) 
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Applying Gauss' theorem to the element and making use of Equations (1.7) and (1.8) gives the 
differential form of Gauss’ theorem: 


1.9 
Ox Oy Ep wea 


OE. OE, OE 
x 4 y 4 so p 
The expression on the left-hand side of Equation (1.9) is known as the divergence of E. It is 
sometimes written as div E. The same expression can also be written as the dot product between the 
differential operator 


cers eee S| 
V=|x—+y—+z (1.10 ) 
Ox “Oy @z 


and the vector E. In Equation (1.10) X, y,Z are unit vectors along the x-, y- and z-axes. Using the 


symbol V , which is known as ‘del’. Equation (1.9) can be written 


ViEs (1.11) 
€ 


This abbreviation is not as pointless as it seems because Equation (1.11) is valid for all systems of 
coordinates in which the coordinate surfaces intersect at right angles. An appropriate form for V can 
be found for each such coordinate system. 


1.5 Electrostatic potential 


The electric field is inconvenient to work with because it is a vector; it would 

be much simpler to be able to work with scalar variables. The electrostatic E 
potential difference (V) between two points in an electric field is defined as 
the work done when unit positive charge is moved from one point to the other. 
Consider the figure on the right. The force on the charge is E, from Equation 
(1.3), so the external force needed to hold it in equilibrium is -E. The work 
done on the charge by the external force when it is moved through a small 
distance dl is the product of the external force and the distance moved in the 
direction of that force. Thus the change in electrostatic potential is 


dV =-E-dl (1.12) 


The potential difference between two points A and B can be calculated by integrating Equation (1.12) 
along any path between them. Mathematically this is written 


B 
V,-V,=—| E-dl (1.13) 
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This kind of integral is called a line integral. This is a slightly tricky concept, but its application is 
limited in practice to cases where the symmetry of the problem makes its evaluation possible. 


The electrostatic potential is analogous to gravitational potential, which is defined as the work done in 
moving a unit mass against gravity from one point to another. The change in the gravitational 
potential depends only upon the relative heights of the starting and finishing points and not on the 
path which is taken between them. We can show that the same is true for the electrostatic potential. 
Figure 1.5 shows a possible path between two points 4 and B in the presence of the electric field due 
to a point charge Q at O. 


0 8+0 


Fig. 1.5 When a charge is moved from A to B in the field of another charge at O the change in 
electrostatic potential is found to be independent of the choice of the path APB. It depends only on the 
positions of the ends of the path. 
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The contribution to the integral of Equation (1.13) from a small movement dl of a unit charge at P is 


ATE yr (1.14) 


But the dot product f-dl is simply a way of writing ‘the component of dl in the radial direction’ 
using mathematical notation, and this quantity is just dr, the change in the distance from O. Thus 
Equation (1.14) can be integrated to give 


V,-V,= [, Sipe c | (1.15) 


> dr = 
4 Argyr ATE) \ 1, Vy 


The potential difference between A and B therefore depends only on their positions and not on the 
path taken between them. By using the principle of superposition we can extend this proof to the field 
of any combination of charges. 


Potential differences are measured in volts. They are familiar to electronic engineers from their role in 
the operation of electronic circuits. It is important to remember that potentials are always relative. 
Any convenient point can be chosen as the zero of potential to which all other voltages are referred. 
Electronic engineers are inclined to speak loosely of the voltage at a point in a circuit when strictly 
they mean the voltage relative to the common rail. It is as well to keep this point in mind. 


It follows from the preceding discussion that the line integral of the electric field around a closed path 
is zero. This is really a formal way of saying that the principle of conservation of energy applies to the 
motion of charged particles in electric fields. In mathematical symbols the line integral around a 
closed path is indicated by adding a circle to the integral sign so that 


f E-dl=0 (1.16) 


The principle of conservation of energy often provides the best way of calculating the velocities of 
charged particles in electric fields. 


1.6 Calculation of potential in simple cases 


In simple cases where the electric field can be calculated by using Gauss’ theorem it is possible to 
calculate the potential by using Equation (1.13). More complicated problems can be solved by using 
the principle of superposition. Since scalar quantities are much easier to add than vectors it is best to 
superimpose the potentials rather than the fields. 
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To show how this is done let us consider Fig 1.6 which shows a cross-sectional view of two long 
straight cylindrical rods each of radius a. The rods are parallel to each other with their centre lines d 
apart. Rod A carries a charge q per unit length uniformly distributed and rod B carries a similar charge 
-q. We wish to find an expression for the electrostatic potential at any point on the plane passing 
through the centre lines of the rods. 


Fig. 1.6 A cross-sectional view of a parallel-wire transmission line. 


The electric field of either rod on its own can be found by applying Gauss’ theorem as described 
above with the result given by Equation (1.6). Since E is everywhere radial it follows that V depends 
only on r and 


f= 4 dy =-—4 _Iny + constant (1.17) 
20E or 2KE, 


where the constant of integration can be given any convenient value because the zero of potential is 
arbitrary. 


It is appropriate to choose the origin of coordinates to be at O, mid-way between the rods, because of 
the symmetry of the problem, so that OA lies along the x-axis. Then for rod 4, r = | x-d / 2 | The 
same argument can be used for rod B, giving an expression for the potential which 1s identical to 
(1.17) except that the sign is reversed and r = | x+d/2 | . Superimposing these two results and 
substituting the appropriate expressions for the radii we get: 

q 


V= In 
278 5 


(1.18) 


wee) 
2x-d 


where the constant of integration has been set equal to zero. This choice makes V = 0 when x = 0. The 
same method could be used to find a general expression for the potential at any point in space. 
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1.7 Calculation of the electric field from the potential 


We have so far been concerned with means of calculating the potential from the electric field. In 
many cases it is necessary to reverse the process and calculate the field from a known potential 
distribution. Figure 1.7 shows how a small movement dl may be expressed in terms of its components 


as 


dl=xdx+ydy+zdz (1.19) 


Fig. 1.7 A small vector dl can be regarded as the sum of vectors dx, dy and dz along the coordinate 
directions 
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The electric field may likewise be expressed in terms of its components 
E=xXE,+yE£,+72£, (1.20) 
Then from Equations (1.12), (1.19) and (1.20) the potential change along dl is given by 


dV = -(KE, +E, +2E,)-(kdx+§dytidz) 
(1.21) 
= - (E,dx + Ey dy +E.dz) 


Setting dy = dz = 0 we have, for a movement in the x-direction, 


with similar expressions for the other two coordinate directions. By superposition the total electric 
field is 


Pane ee ae (1.22) 
Ox Oy Oz 


The expression in parentheses on the right-hand side of Equation (1.22) is termed the gradient of V. It 
can be obtained by operating on V with the operator V defined by Equation (1.10) so that Equation 
(1.22) may be written 


E=-— gradV =-VV (1.23) 


This equation, like Equation (1.11), can be written in terms of any orthogonal coordinate system by 
using the appropriate form for VV . 


From Equation (1.12) it can be seen that if dl lies in such a direction that V is constant it must be 
perpendicular to E. Surfaces on which V is constant are known as equipotential surfaces or just 
equipotentials. They always intersect the lines of E at right angles. It has already been mentioned that 
field plots are useful aids to thought in electrostatics. They can be made even more useful by the 
addition of the equipotentials. Figure 1.8 shows, as an example, the field plot for the parallel wires of 
Fig. 1.6. 
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Fig. 1.8 The field pattern around a parallel-wire transmission line 


1.8 Conducting materials in electrostatic fields 


A conducting material in the present context is one which allows free movement of electric charge 
within it on a time scale which is short compared with that of the problem. Under this definition 
metals are always conductors but some other materials which are insulators on a short time scale may 
allow a redistribution of charge on a longer one. They may be regarded as conducting materials in 
electrostatic problems if we are prepared to wait for long enough for the charges to reach equilibrium. 
The charge distribution tends to equilibrium as exp(- t/ t) , where the time constant t is known as the 


relaxation time. Some typical values are: 


copper 15x10" s 
distilled water 10° s 
fused quartz 10°s 


Once the charges have reached equilibrium there can be no force acting on them and the electric field 
within the material must be zero. 


When an uncharged conducting body is placed in an electric field, the free charges within it must 
redistribute themselves to produce zero net field within the body. Consider, for example, a copper 
sphere placed in a uniform electric field. The copper has within it about 10°’ conduction electrons per 
cubic metre, and their charge is balanced by the equal and opposite charge of the ionic cores fixed in 
the crystal lattice. The available conduction charge is of the order of 10'° C m®, and only a tiny 
fraction of this charge has to be redistributed to cancel any practicable electric field. This 
redistribution gives rise to a surface charge, somewhat as shown in Fig. 1.9, whose field within the 
sphere is exactly equal and opposite to the field into which the sphere has been placed. 
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1. Electrostatics in free space 


—— 


Direction 
of the 
applied 

F field 


Fig. 1.9 The field pattern of the charge induced on a conducting sphere placed in a uniform electric 


field. 
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This surface charge is known as induced charge. It is important to remember that the positive and 
negative charges balance so that the sphere still carries no net charge. The complete solution to the 
problem is obtained by superimposing the original uniform field on that shown in Fig. 1.9 to give the 
field shown in Fig. 1.10. Note that the flux lines must meet the surface of the sphere at right angles 
because the surface is an equipotential. 


ge 


= cia 


Fig. 1.10 The field pattern around a conducting sphere placed in a uniform electric field. This pattern 
is obtained by superimposing the field of the induced charges (shown in Fig. 1.9) on the uniform field. 


Not only is there no electric field within a conducting body, but 

there is also no field within a closed conducting shell placed in 

an electric field. To prove this, consider the figure on the right, 

which shows a closed conducting shell S|. This must be an 

equipotential surface. If there is any electric field within S| there 

must be other equipotentials such as S) lying wholly within 5}. S, 
Now the interior of the shell contains no free charge so, applying 

Gauss’ theorem to 5S, the flux of E out of S> is zero. But, since it 

has been postulated that S> is an equipotential surface, this can be 

true only if E is zero everywhere on it and the potential of S> is the same as that of 5). 


A closed hollow earthed conductor can therefore be used to screen sensitive electronic equipment 
from electrostatic interference. The screening is perfect as long as there are no holes in the enclosure, 
for example to allow wires to pass through. Even when there are holes in the enclosure the screening 
can still be quite effective, for reasons which will be discussed in the next chapter. When the electric 
field varies with time other screening mechanisms come into play and the screening is no longer so 
perfect. 
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1.9 The method of images 


We have already seen that the electric field produced by a known distribution of charges can be 
calculated in simple cases, by the application of Gauss’ theorem and the principle of superposition. In 
most practical problems, however, the charge is unknown and the problem is specified in terms of the 
potentials on electrodes. Simple problems of this type can be solved by the use of Gauss’ theorem if it 
is possible to make assumptions about the distribution of charges from the symmetry of the problem. 


If an uncharged, isolated, conducting sheet is placed in an electric field, then equal positive and 
negative charges are induced on it. Normally this process requires currents to flow in the plane of the 
sheet, and the field pattern is changed so that the sheet becomes an equipotential surface. If, however, 
the sheet is arranged so that it coincides with an equipotential surface, the direction of current flow is 
normal to the plane of the sheet and the two surfaces become oppositely charged. If the sheet is thin, 
the separation of the positive and negative charges is small and the field pattern is not affected by the 
presence of the sheet. This fact can be used to extend the range of problems which can be solved by 
elementary methods. For example, a conducting sheet can be placed along the equipotential AB in Fig. 
1.8. It screens the two charged wires from each other so that either could be removed without 
affecting the field pattern on the other side of the sheet. Thus the field pattern between a charged wire 
and a conducting plane is just half of that of a pair of oppositely charged conducting wires. 
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The field between a charged wire and a conducting plane can be found by reversing the train of 
thought. We note that an image charge can be placed on the opposite side of the plane to produce a 
field which is the mirror image of the original field. The image charge is equal in magnitude to the 
original charge, but has the opposite sign. The plane is an equipotential surface in the field of the two 
charges, so it can be removed without altering the field pattern. The problem is then reduced to the 
superposition of the fields of the original and image charges. This method is known as the method of 
images. It can be applied to the solution of any problem involving charges and conducting surfaces if 
a set of image charges can be found such that the equipotentials in free space of the whole set of 
charges coincide with the conducting boundaries. 


1.10 Laplace’s and Poisson’s equations 


The method described in the previous section has been applied with ingenuity to a wide variety of 
problems whose solutions can be looked up when required. Unfortunately engineers are not free to 
choose the problems they wish to solve, and the great majority of practical problems cannot be solved 
by elementary methods. Figure 1.11 shows a typical problem: an electron gun of the kind used to 
generate the electron beam in a microwave tube for satellite communications. 


Cathod 
eereey Electron 


Focus 

electrode 
at cathode 
potential 


oS Equipotentials 


Fig. 1.11 The arrangement of a typical high-power electron gun. Such a gun might produce a 50 mA 
electron beam 2 mm in diameter for a potential difference between cathode and anode of 5 kV. 


In this case the field problem and the equations of motion of the electrons must be solved 
simultaneously because the space charge of the electrons affects the field solution. A general method 
which can be used, in principle, to solve any problem is obtained by combining Equations (1.9) and 
(1.22) to give 


OV av . av _ Pp 


1.24 
Ox? Oy? zx’ Ey : 
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This is known as Poisson’s equation. It can also be written 


Vres! (1.25) 


a» . : P 
where V~ is given, in rectangular Cartesian coordinates, by 


a on o 
Vs —stas ta (1.26 ) 
Ox” Oy Oz 
When there is no free charge present the equation takes the simpler form known as Laplace’s 
equation: 
2 2 2 
es ga a cae (1.27) 


+ = 
6x? Oy? G2? 


This equation has been solved for a very wide range of boundary conditions by analytical methods 
employing a variety of coordinate systems and by the special method known as conformal mapping, 
which applies to two-dimensional problems. These solutions can be looked up when they are required. 
Cases whose solutions are not available in the literature must, in nearly every case, be solved by 
numerical methods. When free charges are present in a problem it is necessary to use Poisson’s 
equation as the basis of either an analytical or a numerical solution. There are only a few cases which 
can be solved analytically. 


In every kind of active electronic device electric fields are used to control the motion of charged 
particles. The methods described here can be applied to the motion of charged particles in vacuum. 
When the charge densities are small it is possible to calculate the electrostatic fields, neglecting the 
contributions of the charges to them, and then to integrate the equations of motion of the particles. At 
higher charge densities the fields are affected by the space charge and it is necessary to find mutually 
consistent solutions of Poisson‘s equation and the equations of motion. The motion of charge carriers 
in semiconductor devices such as transistors requires knowledge of the fields in material media as 
discussed in Chapter 2. 


1.11 The finite difference method 


The simplest numerical method for solving field problems is the finite difference method. In this 
method a regular rectangular mesh is superimposed upon the problem. The real continuous variation 
of potential with position is then approximated by the values of the potential at the intersections of the 
mesh lines. Figure 1.12 shows a small section of a two-dimensional mesh with a spacing / in each 
direction and the electrostatic potentials at the mesh points. 
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Fig. 1.12 Basis of the finite difference calculation of potential. 
To find an approximate relationship between the potentials shown we apply Gauss' theorem to the 


surface shown by the broken line. The component of the electric field normal to the section AB of the 
surface is given approximately by 


E nap > V, ~ V, yh (1.28) 
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The flux of E through unit depth of the face AB is therefore 


O =: Fok 
AB nAB (1.29) 


If the Gaussian surface does not enclose any charge the net flux of E out of it must be zero so that 


AV, —¥, he —¥ -K, =0 


1 
Vs =i +V,+V,+V,) (1.30) 


Thus, if we know the potentials at points 1 to 4 approximately, we can use (:) 
Equation (1.30) to obtain an estimate of Vo. Because the errors in the four 


potentials cancel each other out to some extent, and because the resulting error G) (+) C) 
is divided by 4, the error in the value of Vp is normally less than the errors in 


the potentials used to calculate it. Equation (1.30) is conveniently summarized 
by the diagram on the right. 


10 mm 


This method can be used to find the fields around two-dimensional 


arrangements of electrodes on which the potentials are specified such as 


the concentric square tubes shown in the figure on the right. The method 16 cur: 


can be implemented on a spreadsheet as follows: 


a) 


b) 


c) 


d) 


A uniform square mesh is defined such that the electrodes coincide 

with mesh lines. The mesh spacing is chosen so that it is small enough to provide a reasonably 
detailed approximation to the fields whilst not being so small that the computational time is very 
large. 


Cells of the spreadsheet are marked out such that one cell corresponds to each mesh point. The 
symmetry of the problem can be used to reduce the number of cells required. Thus, for the 
geometry shown above it is sufficient to find the solution for one quadrant of the problem. 


The electrode potentials are entered into the cells corresponding to the electrodes and the formula 
in Equation (1.30) is entered into all the other cells. When symmetry has been used to reduce the 
size of the problem the formulae in the cells along symmetry boundaries make use of the fact that 
the potentials on either side of the boundary are equal. 


The formulae in the cells are then applied repeatedly ( a process known as iteration ) until the 
numbers in the cells cease to change. To do this the calculation options of the spreadsheet must be 
set to permit iteration. The final numbers in the cells are then approximations to the potentials at 
the corresponding points in space. 
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e) From this solution the equipotential curves can be plotted by interpolation between the potentials 
at the mesh points and the field components can be calculated at any mesh point. 


The method can be applied to more complicated problems including those with curved electrodes 
which do not fit the mesh and three-dimensional problems. Further information can be found in the 
literature. 


1.12 Summary 


In this chapter, starting from the inverse square law of force between two charges, we have derived a 
range of methods for solving practical problems involving electric fields in free space. The concepts 
of electric field, flux density and potential have been shown to be useful for these purposes. The ideas 
contained in this chapter find their direct application in problems about voltage breakdown between 
electrodes in air and those dealing with the motion of charged particles in vacuum. 


The very limited range of problems which can be solved by elementary methods can be extended by 
the use of the principle of superposition and the method of images. In most real problems, however, 
the electric field can be calculated only by solving Laplace’s or Poisson’s equations. Cases which 
have not been solved before generally have to be tackled using numerical methods. 
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2. Dielectric materials and capacitance 


Objectives 


e To discuss how and why an electric field is affected by the presence of dielectric materials 

e To introduce the electric flux density vector D as a useful tool for solving problems involving 
dielectric materials. 

e To derive the boundary conditions which apply at the interface between different dielectric 
materials. 

e To introduce the idea of capacitance as a general phenomenon which is not restricted to 
capacitors. 

e To demonstrate the calculation of capacitance by the use of Gauss’ theorem, field solutions 
and energy methods. 

e To introduce the idea of stored energy density in an electric field. 

e To discuss the causes of electrostatic interference and techniques for reducing it. 
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2.1 Insulating materials in electric fields 


Very many materials do not allow electric charges to move freely within them, or allow such motion 
to occur only very slowly. These materials are not only used to block the flow of electric current but 
also to form the insulating layer between the electrodes of capacitors. In this context they are known 
as dielectric materials. By making an appropriate choice of dielectric material for a capacitor it is 
possible to reduce the size of a capacitor of given capacitance or to increase its working voltage. 


If a dielectric material is subjected to a high enough electric field it becomes a conductor of 
electricity, undergoing what is known as dielectric breakdown. This controls the maximum working 
voltage of capacitors, the maximum power which can be handled by coaxial cables in high-power 
applications such as radio transmitters, and the maximum voltages which can be sustained in 
microcircuits. It is not always appreciated that because dielectric breakdown depends on the electric 
field strength it can occur when low voltages are applied across very thin pieces of dielectric material. 


In order to understand the behaviour of dielectric materials in electric fields it is helpful to make a 
comparison with that of conductors. Figures 2.1 (a) and (b) show respectively a conducting sheet and a 
dielectric sheet placed between parallel electrodes to which a potential difference has been applied. 
The potential difference is associated with equal and opposite charges on the two electrodes. 


(a) (b) 


Fig. 2.1 Comparison between (a) a conductor and (b) an insulator placed in an electric field. 


The conducting sheet of Fig. 2.1 (a) contains electrons which are free to move and set up a surface 
charge which exactly cancels the electric field within the conductor in the manner discussed in 
Chapter 1. The electrons in the dielectric material, on the other hand, are bound to their parent atoms 
and can only be displaced to a limited extent by the applied electric field. This displacement, however, 
is sufficient to produce some surface charge and the dielectric is then said to be polarized. The 
surface charge is not sufficient to cancel the electric field within the sheet, but it does reduce it to 
some extent, as shown in Fig. 2.1 (b). 


Polarization may also produce a volume distribution of charge, but we shall assume that this does not 
occur in the materials in which we are interested. It is important to remember that the surface charge 
produced by the polarization of a dielectric is a bound charge which, unlike the surface charge 
induced on a conductor, cannot be removed. The polarization charges must also be carefully 
distinguished from any free charge which may reside on the surface of a dielectric. 
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On materials which are good insulators, free charges may persist for long periods and strong electric 
fields may build up as a result of them. These phenomena have many important practical 
consequences, but they are not easy to study theoretically because the distribution of the charges is 
usually unknown. In metal oxide semiconductor (MOS) integrated circuits, for example, it is possible 
for charges to build up on the gate electrodes if they are left unconnected. The electric field produced 
by these charges can be strong enough to cause dielectric breakdown of the silicon dioxide layer. This 
is why special precautions have to be taken when handling these circuits. In what follows we shall 
assume that the dielectric is initially uncharged and that any surface charge is the result of 
polarization. 


To put this subject on a quantitative basis, let us suppose that the electrodes in Fig. 2.1 (b) carry a 
surface charge o per unit area and that the surface charge on the dielectric is o, per unit area. Now, 


assuming that the electric field is everywhere uniform and normal to the electrodes, it can be shown 
that the field outside the dielectric is given by 


E,=— (2.1) 


This result is obtained by applying Gauss’ Theorem to the field between plane, parallel, electrodes. 
Similarly the field within the dielectric is 


E,=——— (2.2) 


z, (= |e- 0 oS 


where €, = o/ (c -O x is known as the relative permittivity of the material. Since 0, <o it 


follows that ¢, > 1. It is unfortunate that the symbol ¢, has been adopted for this property of dielectric 
materials because there could be some confusion between it (a dimensionless quantity) and the 
permittivity, defined as ¢ = ¢¢,, and measured in Farads per metre. Care must be taken not to get 
these symbols confused with each other. 


In order to make the theory simpler, we shall assume that o,, is proportional to o and that ¢, is 
therefore a constant. This assumption holds good for many of the materials used in electronic 
engineering, but it is very important to remember that it is not always valid. In particular, for some 
materials, ¢, may depend on: 
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e the strength of the electric field; 

e frequency (if the field is varying with time); 
e the orientation of crystal axes to the field; 

e the previous history of the material. 


Problems involving linear dielectric materials could be solved by calculating the polarization charges 
and finding the fields resulting from both the free charges on the electrodes and the bound 
polarization charges. This would not usually be easy and it is much better to use an approach in which 
the polarization charges are implicit. To do this we introduce a new vector known as the electric flux 
density, which is defined by 


D=cE (2.4) 


In the example given above the electric flux density outside the dielectric is D = o and that within the 
dielectric is likewise o. In other words D depends only on the free charges, unlike E, which depends 
on the polarization charges as well. 


It can be shown that, subject to the validity of the assumption that ¢, is a constant, the argument given 
above can be generalized to cover pieces of dielectric material of any shape. Gauss’ theorem 
(Equation (1.5)) can thus be written in a form which is valid for problems involving dielectric 
materials: 
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ff D-ds= [J] pdv (2.5) 


or, in differential form 
V-D=p= (2.6) 


2.2 Solution of problems involving dielectric materials 


Many problems in electrostatics deal with sets of electrodes together with dielectric materials. When 
the symmetry of a problem is simple it is possible to use Gauss' theorem in much the same way as in 
Chapter 1. 


Figure 2.2 shows a coaxial cable in which the space between the conductors is filled with a dielectric 
material of permittivity ¢. We wish to find an expression for the electric field within the dielectric 
when the potential difference between the electrodes is Vo. 


é 


Fig. 2.2 The arrangement of coaxial cylindrical electrodes, an idealization of a coaxial cable. 


Assume that the conductors carry charges + g per unit length, with the inner conductor being 
positively charged. Applying Gauss’ theorem as given in Equation (2.5) to a cylindrical Gaussian 
surface of radius 7 we find that the radial component of D is given by 


q 
D.=— 2.7 
" nr “) 


then, using Equation (2.4) 


po! (2.8) 
€ 2ner 


The potential difference between the conductors is found by using Equation (1.13) 


b b qd 
V,-V,=-[" Edr =-[’ Lar =-y, (2.9) 
a a Qanr 
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Then 


b 

V,=—4-1n| > (2.10) 
2mé \a 

The charge gq which was assumed for the purposes of the solution need not be calculated explicitly and 

can be eliminated to give the radial component of the electric field as a function of radial position: 


V 
a 2.11 
rIn(b/a) od 


2.3 Boundary conditions 


When two or more dielectric materials are present it is necessary to treat each region separately and 
then to apply the appropriate boundary conditions at the interfaces. There are three of these 
conditions relating to V, E and D respectively. 


The electrostatic potential is continuous at a boundary, that is, it cannot change suddenly there. The 
physical reason for this condition is that an abrupt change in the potential would imply the presence of 
an infinitely strong electric field. 


To find the boundary condition for the electric field we 
consider an infinitesimal closed path as shown in the figure on 
the right. The path is chosen so that it crosses the boundary 


between two dielectrics having permittivities ¢, and ¢) as omar be a 
shown. If the loop is made very thin, then the contributions to ae | de es 
the line integral in Equation (1.16) arising from the parts of the 

loop normal to the boundary are negligible. If, in addition, the 

tangential components of the electric field are E,; and Ey, then the integral becomes 


(Ey ~E,,)dl =0 
so that 
Ey, =£, (2.12) 


This result can be stated in words as: the tangential component of Eis continuous at a boundary. 


The boundary condition for the electric flux density can be 
found in a similar way by using Gauss’ theorem. The figure in 
the margin shows a boundary between two dielectric materials 
with an infinitesimal Gaussian surface which crosses it. If the 


38 


Electromagnetism for Electronic Engineers 2. Dielectric materials and capacitance 


thickness of the ‘pill box’ is very small, then the contributions to the flux from the parts of the surface 
which are normal to the boundary are negligible and the integral becomes 


(D,, ~ Dy» )dA = 0 


where D,,,, and D,,. are the components of D normal to the boundary on each side of it and dA is the 
area of the part of the boundary lying within the Gaussian surface. Thus 


D = D,y (2.13) 


or, in words, the normal component of D is continuous at a boundary. 

To solve problems with two, or more, layers of dielectric material we first apply Gauss’ theorem to 
find D everywhere since that does not depend on the presence of the materials. Next Equation (2.4) is 
used to find E in each region. Finally the potential difference across each layer is found using 
Equation (1.13). 
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2.4 Capacitance 


Capacitors are very familiar as circuit elements, but it is not always realized that the idea of 
capacitance is more general. Capacitance exists between any pair of conductors which are electrically 
insulated from each other. Thus there is a capacitance between adjacent tracks on a printed circuit 
board, but this does not usually appear in the circuit diagram. This ‘stray’ or ‘parasitic’ capacitance 
can cause unwanted coupling between the parts of a circuit, causing it to oscillate or misbehave in 
some other way. Very few electronic engineers ever need to calculate the capacitance of a capacitor; 
they are much more likely to need to estimate the magnitude of a stray capacitance. 


The figure on the right shows a cross-sectional view of two 

adjacent tracks on a printed circuit board. Let us suppose that the ate -a , 
tracks are insulated from each other and from earth and that S/S TZ/LLD 
charge Q is transferred from B to A. Electrode B must then carry 

charge -Q. As a result of this transfer of charge, an electric field 

exists around the electrodes such that the potential difference between them is V. 


If the dielectric material of the printed circuit board has a permittivity which does not vary with the 
electric field, then the system is linear and the principle of super-position may be applied. It follows 
that the potential difference between the electrodes is directly proportional to the charges on them, so 
we can write 


Q=CV (2.14) 


where C is a constant of proportionality which is readily recognized as the capacitance between the 
electrodes familiar from elementary circuit theory. The unit of capacitance is the farad (F) and 

1 F=1(C V". Most capacitances are small and measured in microfarads, nanofarads or picofarads. 
Stray capacitances are usually of the order of picofarads. 


2.5 Electrostatic screening 


It has already been noted that unwanted capacitive coupling between electronic circuits can be a major 
problem. This is part of the larger problem of electromagnetic interference; another aspect, inductive 
coupling, is discussed in Chapter 6. The problem with all types of electromagnetic interference is how 
to minimize it rather than how to calculate its magnitude accurately. Electromagnetic theory provides 
the means for understanding the causes of electromagnetic interference and the techniques for dealing 
with them. 


A simple case of the coupling of two circuits by stray capacitance is shown in Fig. 2.3. The circuits | 
and 2 are linked by the stray capacitance C, and by a common earth. The stray capacitance is small, 


typically of the order of a pico-farad, so its impedance (Z, = 1/j@C, ) is high, but decreases with 


increasing frequency. In this problem V is the source of the interference picked up by circuit 2. 
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V, 


Fig. 2.3 Showing how two circuits can be coupled together by a stray capacitance C;. The unwanted 
signal coupled from circuit 1 into circuit 2 can be large enough to be troublesome. 


The current flowing through the capacitor is small compared with that in R,) so the spurious signal 
appearing at the input of the amplifier is approximately 


R R' 
i= _ mV, (2.15) 
Ry, +r Ri Zz oF R,, 


where R;, = (R,R,, y (R,, ase i The spurious voltage given by Equation (2.15) is to be compared 


In 


with the signal voltage at the input of the amplifier 


V SRR, +R, ) 


sig in 


Equation (2.15) shows that the spurious signal is greatest when the source impedance Rg), is low and 
the effective input impedance of the second circuit R',, is high. An indication of the order of 
magnitude of capacitance which can cause trouble can be obtained by supposing that circuit | is the 
a.c. main. R,; is then very small. If R'i, 1s 1 MQ and V, is 1 wV, Z, is approximately 2.4 x 10!" Q, 
which corresponds at 50 Hz to a stray capacitance of the order of 10°" F. 


The capacitive coupling between the circuits can be reduced by putting an earthed screen between 
them, as shown in Fig. 2.4. The stray capacitance is divided into two parts in series with each other 
and with their common point earthed. In practice, unless the screen completely encloses one of the 
circuits, there is still a residual capacitance connected directly between P and Q, bypassing the screen. 
Quite effective screening can be achieved with a partial enclosure provided that it intercepts most of 
the field lines passing from P to Q. It is essential that the screen is earthed, otherwise P is connected 
to O by Cy) in series with C,, without a signal path to earth from their common point. 
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Screen Pat 


Fig. 2.4 The unwanted coupling between the two circuits illustrated in Fig. 2.3 can be greatly reduced 
by putting an earthed screen between them. Any current passing through Cs, is conducted to earth 
instead of passing through Cs to O. 


When very good screening against low-frequency electric fields is required then a closed conducting 
box must be used. Since all the sources of the field lie outside the box and it must be an equipotential 
surface it follows that, theoretically, the electric field inside it is zero. Such a box is sometimes known 
as a ‘Faraday cage’. In practice, any enclosure will have one or more holes in it to allow wires to pass 
in and out and these may reduce the effectiveness of the screen. Similarly, the screening effectiveness 
is affected by the way in which the joints of the box are made. 
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The electric screening effectiveness of an enclosure is defined as the ratio of the magnitude of the 
electric field with the screen to that at the same point when the screen is removed. It is usually 
expressed in decibels. At high frequencies other factors come into play. 


The effects of capacitive coupling can also be reduced by using a differential amplifier if the source 
V2can be isolated from earth. Figure 2.5 shows how this works. The stray capacitances C,; and Cs» are 
often approximately equal and so the spurious signals appearing at the normal and inverting inputs of 
the amplifier are also very nearly the same. If the amplifier has a high common-mode rejection ratio, 
then only their difference is amplified and added to the wanted signal. Adding an earthed screen 
reduces the stray capacitive coupling and therefore reduces the unwanted signal still further. The 
source for circuit 2 may be a transducer situated some distance from the amplifier. In that case the 
connecting cables must be screened as well as the amplifier, a point which is discussed further in 
Chapter 6. 


Fig. 2.5 When the input to a circuit can be isolated from earth it is possible to reduce capacitive 
coupling effects by using an amplifier with a differential input. The unwanted signals coupled through 
Cs, and C,p are approximately equal but only their difference is added to the signal. 


It is also possible for there to be capacitive coupling between the output and the input of an amplifier. 


When that happens there can be positive feedback causing the circuit to oscillate, so once again 
screening is needed. 


2.6 Calculation of capacitance 


In simple cases where the field problem can be solved by the direct application of Gauss’ theorem it is 
straightforward to calculate the capacitance. The sequence of steps is: 


1. Assume charges + q on the conductors. 

2. Apply Gauss’ theorem, superposition or the method of images to find E. 

3. Integrate E along any convenient path between the conductors to find the potential difference 
between them. 

4. Calculate the capacitance from C = q/V. 


When the shapes of the electrodes are more complicated than those in the examples above it is no 
longer possible to use the same method to calculate the capacitance. An alternative method is: 
1. Assume potentials 0 and V on the electrodes. 
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2. Solve Laplace’s equation analytically or numerically to give values of the potential 
everywhere. 

3. Calculate E close to the surface of one of the electrodes using E = -grad V. 
Calculate the surface charge density distribution on the electrode using o = éE. 

5. Integrate the charge density distribution over the surface of the electrode to give the total 
charge q. 

6. Calculate the capacitance from C = q/V. 


To summarize: capacitance can be obtained by calculating either 


a) The potential difference between the electrodes from assumed charges, or 
b) The charges on the electrodes from an assumed potential difference between them. 


By using both approaches in an approximate way it is possible to obtain upper and lower bounds to 
the capacitance. The formal application of this in the energy method is described later in the chapter. 


2.7 Energy storage in the electric field 


When charge is transferred from one plate of a capacitor to the other work is done against the electric 
field. This work is stored as potential energy in the capacitor. It is easy to show that the stored energy 
is given by 


2 
Mept © oy (2.16) 


W = = = 
2 2C 2 


It is important to remember that Equation (2.16) applies to any pair of electrodes with a potential 
difference between them, not just to capacitors as lumped components. 


Fig. 2.6 The electric field between two charged conductors can be thought of as being made up of 
tubes of flux like the one shown. The walls of the tube are everywhere parallel to the local direction of 
the electric field. Its ends are terminated by equal and opposite charges. 


It is sometimes useful to think of the energy stored in a capacitor as being distributed throughout the 
electric field associated with it. Figure 2.6 shows a capacitor with electrodes of arbitrary shape. We 
define a flux tube by considering a small element of the surface of the conductor 4 and following the 
lines of E which start from its boundary through space until they terminate upon B. Since the flux 
lines can never meet or cross each other it follows that the whole electric field of the capacitor can be 
divided up into flux tubes. In addition we observe that, if the charge on the element of A from which 
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the tube starts is +dQ, then that on the element of B on which it ends must be -dQ by a straightforward 


application of Gauss’ theorem to a Gaussian surface enclosing the tube. 


The figure on the right shows a small part of the flux tube on a larger 
scale. A short length of the tube is chosen having length d/ as shown, 
and the cross-sectional area of the tube at this point is dA. If the local 
electric field strength is E, then the potential difference between the 
ends of this element of volume is dV = E dl, where the dot product can 
be omitted because E and d/ are parallel to each other by definition. 
Likewise, the relationship between the electric flux in the tube and the 
charges at its ends means that dO = D dA. 


The energy stored in the element of volume is then 


dW = ;doav = SDE dldA = + DEdv 


using Equation (2.16) with the volume of the element represented by dv. In the limit as dv shrinks to 


zero the energy density in the electric field is ’2DE. A more rigorous derivation of this result which 


allowed for the possibility that the vectors D and E are not parallel to each other would show that the 


energy density in the field is given by 
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w=—D-E 2.17 
5 (2.17) 


and the total energy stored in the field of the electrodes is 
1 
W=—|\[ DE dv (2.18) 


where the integral is taken over all that part of space in which the electric field is not zero. 


2.8 Calculation of capacitance by energy methods 


In many cases the solution required to a field problem is not detailed field information but a single 
number, the capacitance. The use of the concept of stored energy provides a particularly simple way 
to obtain estimates of capacitance without finding the exact solution to the field problem. The starting 
point is the fact that when a physical system is in stable equilibrium its stored energy is normally a 
minimum. From this it follows that any perturbation of the system must result in an increase in the 
stored energy. 


Consider now two possible perturbations of the system of charges and electric fields in a capacitor. In 
the first case we retain the electrodes as equipotential surfaces but alter the charges on them in such a 
way that the equipotentials in the inter-electrode space assume particularly simple shapes. From these 
perturbed equipotentials we calculate the energy density in the field and, hence, the stored energy W’ 
of the system. This energy must be greater than the stored energy W of the unperturbed system, so 
that 


(2.19) 


This provides an upper bound to the capacitance. 


In the second case we hold the charges on the electrodes constant but redistribute them in such a way 
that the flux lines have a particularly simple form. When this is done the electrodes are no longer 
equipotential surfaces. Once again the stored energy, say W", is calculated, and this time 


2 2 
w"> Qo or Ce Q (2.20) 


giving a lower bound for the capacitance. 
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2.9 Finite element method 

Many computer packages for the solution of field problems are based upon the use of finite elements. 
This technique is not suitable for hand calculation and the algorithm is a little complicated so only the 
principles of the method as applied to two-dimensional problems are described here. 

In the finite element method the region between the electrodes is divided into a large number of 
triangular elements. Figure 2.7(a) shows one such element in the x-y plane. It can be shown that the 
condition that Laplace’s equation should be obeyed inside the element is equivalent to the 


requirement that the energy stored within it should be a minimum. If the potential within the element 
is approximated by 


V(x,y)=a+bx +cy (2.21) 
then the electric field components are 
E.=-b and E,=-c (2.22) 


and the energy density is 
_l ( 2 2 
we bo +e (2.23) 


Now the coefficients b and c can be expressed in terms of the positions of the vertices of the triangle 
and of the potentials at them. Thus the stored energy in the element can be expressed in terms of the 
same quantities. 


Fig. 2.7 The finite element method: (a) a single element and (b) a cluster of elements. 


Now consider the cluster of triangular elements with a common vertex at A, as shown in Fig. 2.12(b). 
The energy stored in the cluster of elements, W, is the sum of the energies stored in each. We now 
select the value of the potential at A which minimises the stored energy by requiring that 
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i =0 (2.24) 
OV, 

The result is an equation giving the potential at A in terms of the potentials at the other vertices and 
their positions. It is interesting that, if triangular elements are formed from the square elements in Fig. 
1.12 by adding diagonals from the lower left to the upper right corner of each square, then the 
application of the finite element method yields Equation (1.30). 


The finite element method is important chiefly because of the freedom with which the sizes and 
shapes of the elements may be chosen. It is, therefore, possible to select a set of elements which fits 
the boundaries of the problem with good accuracy and which provides a concentration of small 
elements in regions where the potential is changing rapidly. In two-dimensional problems the spatial 
variation of the potential can be represented as a smooth surface. The finite element method 
approximates this surface by one with many small flat facets rather like the surface of a cut gem. 
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2.10 Boundary element method 


The finite difference and finite element methods both depend upon the division of the space between 
the electrodes into a large number of small elements. The potential values at the nodes of the mesh are 
manipulated according to some algorithm until the solution has converged. The result of the process, 
in each case, is detailed information about the field at every point within the solution space. In large 
three-dimensional problems the number of variables can be very large indeed and the computation 
time correspondingly long. 


Very often the result required from a field calculation is not the detailed field information but only a 
single number such as the capacitance. For problems of this kind an alternative method, the boundary 
element method, is useful because it only requires mesh points and potentials to be specified on the 
boundaries. This has the benefit that the number of variables in the solution is less than in the 
corresponding finite difference or finite element solution by two or three orders of magnitude. It also 
makes it easier to model problems with complicated boundaries. The fields at internal points can be 
calculated from the solution at the boundary if required. 


2.11 Summary 


Insulating materials are widely used in electronic engineering both to provide electrical isolation and 
to enhance the performance of capacitors. The dielectric qualities of most insulating materials can be 
represented by a constant, the permittivity. It is convenient to introduce a new field vector, the electric 
flux density D, whose value is independent of the presence of dielectric materials. The use of this 
vector makes the solution of problems involving dielectric materials easier. It is important to use the 
correct boundary conditions at the interfaces between different dielectric materials. 


Capacitance occurs in capacitors and in stray capacitance between conductors. It can be calculated 
either by a direct solution of the field problem or by considering the stored energy in an 
approximation to the field associated with the capacitor. The latter method gives quite good accuracy 
with relative ease. Stray capacitances lead to unwanted coupling between electronic circuits. 
Techniques for reducing the coupling can be understood using the theory of electrostatics. 
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3. Steady electric currents 


Objectives 

e To introduce the ideas of current flow in electrical conductors and of electrical resistance. 

e To show that Ohm’s law is a special case of a more general relationship between current 
density and electric field. 

e To derive expressions for the power dissipated per unit volume within a conductor. 

e To show how the distribution of current within a conductor can be calculated. 

e To introduce the continuity equation and to show how Kirchhoff’s current law is a special 
case of it. 

e To introduce the concept of the electromotive force of a source of electric power. 

e To modify the equation for the line integral of the electric field around a closed path derived 
in Chapter | to allow for the presence of sources of electric power, and to show that the 
resulting equation is a generalization of Kirchhoff's voltage law. 

e To show how the resistance of a conductor can be calculated directly, by the use of Laplace’s 
equation, and by energy methods. 


3.1 Conduction of electricity 


In the preceding chapters we have studied problems involving charges which are at rest. This is 
convenient for the formal exposition of electromagnetic theory, but most people meet electrical 
phenomena for the first time in the form of electric currents flowing in wires. A current can flow in 
any medium in which there are charges which are free to move. These conduction charges may be 
electrons, positively charged ‘holes’, or positive or negative ions, according to the material. They are 
in continuous random motion, with a distribution of velocities which depends upon the temperature of 
the material. They are also constantly colliding with each other and with the atomic structure of the 
material. 


When a conducting material is placed in an electric field the conduction charges are accelerated in the 
direction of the field. The velocity acquired is small compared with the average value of the random 
velocity at ordinary temperatures, and is superimposed on it. The ordered part of the motion would 
increase without limit were it not for the collisions which convert it into random motion. Overall, the 
effect of the field is to add a small average drift velocity to the random velocity. The magnitude of the 
drift velocity is related to the strength of the field by the equation 


v,=HE (3.1) 


where yz is known as the mobility of the charge carriers. The mobility depends on 
1. the type of charge carrier; 
2. the material in which it moves; 
3. the temperature of the material. 

It can also depend on 


50 


Electromagnetism for Electronic Engineers 3. Steady electric currents 


4. the strength of the electric field; 

5. the orientation of the crystal axes to the field. 
For many materials is a scalar quantity which varies only with temperature and this will be assumed 
here. 


The rate of flow of electric charge across unit area of a plane normal to the direction of vq is the 
current density given by 


J=nqv, (3.2) 


where v1 is the density of charge carriers and q the charge on each. Combining Equations (3.1) and 
(3.2) gives 


J=nquE=cE or E=p—J (3.3) 


where o is the conductivity of the material, measured in siemens per metre. J is measured in 
coulombs per square metre per second or in amperes per square metre, defining the ampere as a 
current of one coulomb per second. It is sometimes convenient to use the reciprocal of conductivity 
which is known as resistivity (). Resistivity is measured in V m A’ or Q m. 
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To link these ideas with those of elementary circuit theory Cross-sectional 
consider the flow of electric current in a section of straight ahead 
wire as shown in the figure on the right. If the wire has cross- — 
sectional area A and carries a current / uniformly distributed , 

across it, then the current density is hoe | 


Vv 


J=1/A (3.4) 


The relationship between the electric field in the wire, the potential difference between its ends, and 
its length is 


E=V/l (3.5) 


But EF = pJ from Equation (3.3) and therefore 
pl 

Vs|2oyeRi (3.6) 
A 


which is the familiar form of Ohm’s law. Equation (3.3) is thus shown to be a general form of Ohm’s 
law applicable to problems in which the current density is not uniform. 


3.2 Ohmic heating 


When conduction charges are accelerated by an electric field, they continually gain energy from it. 
This ordered kinetic energy is transferred to the bulk of the material by collisions, so increasing the 
random thermal motions of the atoms. 


The conversion of electrical energy into heat is familiar from its everyday use in electric heaters and 
light bulbs. For a steady current / sustained by an applied voltage V the power input is 


P=VI=D’R=V7[R (3.7) 


an equation which is well known from elementary circuit theory. By using the results of the previous 
section we can find corresponding formulae for the details of the energy dissipation within a 
conducting material. 

& 
Consider the element of volume shown in the figure on the right. Provided a 
that the properties of the material are the same in all directions, we can set 
the x-axis parallel to the direction of the current without loss of generality. 
The potential difference across the element is 
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V=E& (3.8) 
and the current flowing though it is 

Il=J ov& (3.9) 
Thus the power dissipated in the element is 

OP = EJ & by & = EJiv = 0 E’ bv = pJ’ &v (3.10) 
where Ov is the volume of the element. A more rigorous derivation would show that 

OP=E-J ov (3.11) 
This expression is valid even when the vectors E and J are not parallel to each other. Equation (3.10) 
shows that the power density is greatest in regions of high current density. This phenomenon is of 
great practical importance. In fuses the provision of a section of thin wire in a circuit ensures that that 
section is heated to melting point before any other part of the circuit is damaged. On the other hand, 


wiring joints which are carelessly designed or made can overheat and fail. 


3.3 The distribution of current density in conductors 


We would like to know how to calculate the variation of current A 

density with position. To do this we consider the net current s * 

flow out of a closed surface S enclosing a volume Vas shownin wwH_. we 

the figure on the right. If the current density in the regionofSis — __e 

J (not necessarily a constant) then the net current flow out of S 

is given by the integral of the normal component of J over S, ee ae 

that is 

current out of S = ff J-dA (3.12) 
s 


Now if there is a net current flow out of a closed surface, then the total charge enclosed by that 
surface must be changing with time because charges cannot be created or destroyed. Thus 


current out of S = -£ (charge enclosed by S) 


which may be written in mathematical notation as 
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ff s-aa =-£/ ff pdv (3.13) 
S V 


This is a mathematical statement of the principle of conservation of charge. Note that p is charge 


density, not resistivity, in this equation. 


Using the same method as we used in Chapter | to get from Equation (1.5) to Equation (1.11), it is 
easy to show that the differential form of Equation (3.13) is 


vy.ju- (3.14) 
Ot 


This equation, in either its integral or its differential form, is known as the continuity equation. It is 
valid for the flow of charges in a vacuum as well as in conducting materials. 


When the current flow is in a steady state, the right-hand sides of Equations (3.13) and (3.14) are zero. 
If, furthermore, the conductivity of the material is constant, then 


V-J=V-(cE)=0(V-E)=0 (3.15) 
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since o is a constant. But 

E=-VV (Equation (1.23)) 
so that, from Equation (3.15) 
V-(VV)=VV =0 


Thus when a steady electric current flows in a material of constant conductivity the potential 
distribution obeys Laplace’s equation (Equation 1.27). 


In the special case of currents flowing in wires, as shown in Fig. 3.1, the integral on the left-hand side 
of Equation (3.13) becomes the sum of the currents in the wires 


tities. 41, =0 (3.16) 


which is Kirchhoff’s current law familiar from electrical circuit theory. The continuity equation is 
therefore a generalization of Kirchhoff’s current law. 


a 


Fig.3.1 Applying the continuity equation to the special case of steady currents flowing in wires 


In general, provided that the conductivity is constant, the current 

density distribution can be calculated from the solution of Laplace’s 

equation. The figure on the right shows a typical problem: a strip with 

a step change in its width. In this case Laplace’s equation would have 

to be solved by numerical methods. The potential is specified by 

assuming that at some distance from the step the equipotentials are 

straight lines perpendicular to the length of the strip. Along the sides there is a boundary between a 
conductor and an insulator. Since no current can cross this boundary the electric field must be 
tangential to the boundary. 
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3.4 Electric fields in the presence of currents 


When electric currents flow through a system of conductors there is an electric field in the space 
around the conductors in addition to the field within the conductors which drives the current. As an 
example, consider the field around a coaxial line carrying equal and opposite currents in the core and 
the sheath, as shown in Fig. 3.2. Within the conductors the field lines, like the current flow lines, lie 
parallel to the axis. Outside them the field is modified from the purely radial electrostatic solution. 


The finite resistance of the conductors requires that the potential at A should be higher than that at B 
and that at C should, likewise, be greater than that at D. There must, therefore, be axial components of 
the electric field in the space between the conductors as well as in the conductors themselves. The 


tangential component of the electric field is continuous at the boundary (Equation (2.12)) so the field 
lines are curved as shown and are no longer normal to the boundary. 


A A Le a a a a a a, 


} EEEEEEEL 
A Ps 


Fig. 3.2 Electric field lines and equipotential surfaces in a coaxial cable with resistive conductors. 


3.5 Electromotive force 


In this chapter we have assumed, so far, that sources of electric current exist, but we have paid no 
further attention to them. These sources have been represented as idealized voltage generators. In 
practice, they are most likely to be electrochemical cells or electromagnetic generators. The details of 
the inner workings of these sources do not fall within the scope of this book, but it is important to 
consider the effects of sources in the theory of electromagnetism. 


We shall regard a source as a device capable of maintaining a potential difference and a steady flow 
of electric current between its terminals. Such a device is a source of electric power because it does 
work in driving the electric current around the circuit against the circuit resistance. As the charges 
move around the circuit from one terminal of the source to the other, they lose potential energy as it is 
converted into heat. 


In Chapter | we saw that in an electrostatic field the line integral of the electric field around a closed 
path is zero 
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f E-dl=0 (Equation (1.16)) 


This equation remains true near a current-carrying circuit provided that the path of integration does 
not pass through the generator. But if the line integral of the electric field is evaluated around the 
circuit from one terminal of the source to the other, the result is not zero. This is because a net amount 
of work has to be done to move the charges along this path against the resistance of the circuit. The 
potential energy given up by a charge in this way is restored when it moves through the source. This 
can be expressed mathematically by writing 


[ E-dl=~ (3.17) 
CG 


where !! is the electromotive force (e.m.f.) of the source and the integral is taken around the circuit 
(C) from one terminal of the source to the other. 


When the source is imperfect it is usual to equate the electromotive force with the Thévenin voltage, 
i.e. the open-circuit voltage of the source, and to include the source impedance in the external circuit. 
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For the special case of a circuit consisting of lumped components and generators joined by lossless 

conductors, the integral in Equation (3.17) becomes the sum of the potential differences across the 

components. The electromotive force is, likewise, the sum of the e.m.f.’s of all the generators in the 
circuit. The result is Kirchhoff’s voltage law: 


Ri tial, feat = a a (3.18) 


3.6 Calculation of resistance 


From Ohm’s law in the form 
pl ; 
V= 7 IT=RI (Equation (3.6)) 


it is clear that resistances can be calculated either by assuming a current and working out the potential 
difference, or by the opposite process. It is also possible to make use of energy methods, with the upper 
and lower bounds for the resistance being given by these two approaches. In simple cases the resistance 
of a conductor can be calculated directly using the expression in the brackets in Equation (3.6). 


This approach is only of use in the relatively small number of cases where the shape of the conductor 
is such that either the current, or the potential, distribution is simple. This is clearly not the case with 
shapes such as that shown in Fig. 3.3. Problems of this kind can be solved by the following procedure: 

1. Assume a value for the potential difference V between the ends of the strip; 

2. Solve Laplace’s equation with the appropriate boundary conditions; 

3. Calculate the electric field and current density distributions on any convenient line across the 

strip using EF =—VV and J=oE£; 
4. Integrate J across the strip to find the total current J; 
5. R=DV/I. 


These steps can be carried out by numerical methods if necessary. 


Resistivity p 


Fig. 3.3 Electric field lines (shown dashed) and equipotential surfaces for a current flowing through a 
conducting strip with a step change in width. 


This method can sometimes be a bit tricky to apply because of the need to ensure that the normal 
component of E is zero at the edge of the strip. We can sometimes get round the difficulty by using 
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the principle of duality. In a field map such as Fig. 3.3 the field lines and the equipotentials are 
always at right angles to each other. We can therefore imagine a second problem whose field pattern 
is identical to the first except that the roles of the field lines and the equipotentials have been 
exchanged. Two problems related to each other in this way are called duals of each other. Sometimes 
it is easier to solve the dual of a problem and then to deduce the required answer from that solution. 


3.7 Calculation of resistance by energy methods 


It has already been observed that capacitances can be calculated by energy methods to an accuracy 
which is adequate for many purposes. Similar methods can be used to calculate resistances. In this 
case it is assumed that the effect of perturbing the distribution of the electric field and the current 
density from equilibrium is to increase the rate of energy dissipation. Now the power dissipated in a 
resistance can be written as 


W=I?R=V?/R (Equation (3.7)) 


so, if the total current is given, and the power dissipated is calculated as W’ from an approximate 
distribution of the current, then 


W'>IPR sothat R<W'/I? (3.19) 


giving an upper bound for the resistance. If, on the other hand, the potential difference across the 
resistance is given and the power W" is calculated from an approximate set of equipotentials then 


W">V?/R sothat R>V?/W" (3.20) 
giving a lower bound for R. 
3.8 Summary 


In this chapter we have explored the links between the theory of electromagnetism and the theory of 
electric circuits. It has been shown that Ohm’s law and Kirchhoff’s laws are special cases which arise 
when electromagnetic theory is applied to lumped components and to electrical networks. 


On occasion it is necessary to calculate the resistance of a piece of conducting material. Three 
possible approaches have been discussed: direct calculation, the use of Laplace’s equation and the use 
of energy methods. Direct calculation is possible only in a limited number of cases where the current 
density is uniform over surfaces of simple shape. The application of Laplace’s equation can, in 
principle, solve all problems involving uniform conducting materials. Similar numerical methods can 
be used when the materials are not uniform. These methods are available as computer packages of 
great sophistication. For many purposes great accuracy is not required and then the energy methods 
are of value because of their ability to provide acceptable accuracy with little effort. 
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4. The magnetic effects of electric currents 


Objectives 
e To explain the relationship between magnetism and electricity. 
e To introduce the Biot-Savart law and the magnetic circuit law and to show how they are used 
to calculate the magnetic flux density produced by electric currents. 
e To introduce the concept of magnetic scalar potential and to discuss the types of problem in 
which it can be used. 
e To show that forces are exerted on current-carrying conductors in magnetic fields. 


4.1 The law of force between two moving charges 


The theory of magnetism was developed independently from that of electricity until Oersted (1820) 
showed that the two subjects were linked. This link was later explored by Ampere, Faraday and 
Maxwell, but it was not until the theory of relativity was developed by Einstein that the relationship 
between them could be fully understood. It is now clear that magnetic effects can be regarded as a 
consequence of the motion of electric charges. It is not necessary to postulate the existence of 
magnetic poles or dipoles to explain experimental observations. The magnetic properties of materials, 
discussed in Chapter 5, can be explained by assuming the existence of circulating currents on the 
atomic scale. In this chapter we shall confine our attention to magnetic effects in free space. 
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Q 
The figure on the right shows two charges which have - V, 
constant velocities v; and v>. It can be shown that the force P 
exerted on Q> by Q,, is given by V, 
Q, 
par 2 HE wy (vy, at))=E,+E, (4.1) 
Are r Arr 


where F is a unit vector pointing from | to 2 and the symbol A represents the vector cross product. 
The first term in this equation is the electrostatic force (Equation (1.1)). The second term, which we 


shall recognize as the magnetic force, is a result of the motion of the particles. The constant //, is 


known as the primary magnetic constant; it has the value 42 10’ H m". (The Henry (H) is the unit 
of inductance — see Chapter 6). 


It is instructive to compare the magnitudes of the two forces F, and Fi». To simplify matters let us 
suppose that v; is parallel to v2. This gets rid of the vector products without affecting the order of 
magnitude of the result. The ratio of the magnitude of the magnetic force to that of the electrostatic 
force is then 


i VV, 
—“~=¢ UWv.v, =——_+ (4.2) 
F pee Cc 


m 


It can be shown that the constant c is the velocity of light (see Chapter 7), so c’ has a numerical value 
which is approximately 10'’ m’s”. Thus the magnetic force is negligible compared with the 
electrostatic force unless either: 

1. both particles have velocities close to that of light, or 

2. the electrostatic component of the force is cancelled by the electrostatic force of an equal and 

opposite stationary charge. 

The latter is just what happens when electrons are moving along conducting wires because those 
contain equal numbers of fixed positive charges. Although the electrons move at velocities much less 
than that of light, there are enormous numbers of them, so useful forces are produced. That is why 
magnetic forces can be used to generate electric power and to drive electric motors. Note that the 
cross product in Equation (4.1) means that the magnetic force on Q> is always at right angles to vo. 


4.2 Magnetic flux density 


In the development of the theory of electrostatics in Chapter | the sources of electrostatic force were 
separated from their effects by introducing the concept of the electric field. Following the same 
approach, we introduce a new vector B and split the magnetic term F,, in Equation (4.1) into two 
parts: 
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_ MQ, A 
B= real (v, AP) (4.3) 
and 
F,, = Q,(v, AB) (4.4) 


The vector B is referred to, somewhat confusingly, as the magnetic flux density. The units of B are 
tesla (T) or webers per square metre (Wb m”). Equation (4.3) provides the basis for calculating the 
magnetic flux density produced by any combination of charges moving in free space. The effect of 

that flux density on any other charges can then be found by using Equation (4.4). 


For the moment we will concentrate on ways of computing magnetic flux densities. This information 
might be needed to calculate: 

e the forces on moving charges; 

e the forces on current-carrying conductors; 

e the forces on iron surfaces; 

e — self- and mutual inductances; 

e induced voltages and currents; 

e eddy current losses. 
The second and third of these are mainly applied to electric machines, and we shall refer to them only 
briefly. The last two are dealt with in chapter 6. 


The source of the magnetic flux is often a current flowing in 

a wire, so another form of Equation (4.3) is needed. 

Consider an element of length d/ in a wire carrying a current 

Tas shown in the figure on the right. Let the current consist 

of charges of magnitude g moving with mean drift velocity at a 
va. If there are n charge carriers per unit volume and the were “A 
cross-sectional area of the wire is A, then the current in the 

wire is given by J = ngAv, from Equation (3.2). But the 


charge contained in the element of wire isnqA dl, so that 
Ov, =(ngAdl)v, =1dl (4.5) 


The contribution to the magnetic flux density from the charges within the element is therefore, from 
Equation (4.3), 


dB = uolng dl) (y at)= (alae) (4.6) 
Azar 4ar 


where dl is a vector element of length. 


62 


Electromagnetism for Electronic Engineers 4. The magnetic effects of electric currents 


But isolated current elements cannot exist so Equation (4.6) must be integrated around the complete 
circuit in which the current is flowing. The result is 


p= “lg dlar (47) 


4a r 
This equation is known as the Biot-Savart law. It can be used to compute the magnetic flux density at 
any point produced by known currents flowing in a given arrangement of conductors. Since the flux 
density is proportional to the current such problems are linear and the principle of superposition can 
be used. 


Figure 4.1 shows a short straight section of wire carrying a current /. To find an expression for the 


magnetic flux density at the point P at a perpendicular distance R from the wire we consider a small 
element of the wire of length d/. 
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i-~—_- dé 


P 
Fig. 4.1 Calculation of the magnetic flux density produced at the point P by a short length of wire 


carrying current /. 


From Equation (4.6) the magnetic flux density at P produced by the element d/ is 


GB = nae 
47R* sec’ 0 


(4.8) 
But 
dl = d(Rtan@)= Rsec’ 0d 


so, substituting for d/ in Equation (4.8), 


dB = H,/cosA0dé 
4zR 


The magnitude of the flux density at P produced by the whole wire is found by integrating over the 
length of the wire 


2 Isi 
p= | cosod9 = He SY 


(4.9) 
47R -- 27R 


Examination of the vector cross product in Equation (4.6) shows that B is perpendicular to both R and 
to the wire, and that it is directed into the paper. 
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4.3 The magnetic circuit law 


A useful method for calculating the magnetic flux density can be derived by considering the line 
integral of the vector B around a closed path. Figure 4.2 shows such a path encircling a long straight 
wire carrying current /. The return wire for the current is assumed to be far enough away for it to have 
no appreciable field anywhere on the path. The flux density at any point on the path can be calculated 
from the Biot-Savart law by allowing a to tend to 2/2 in Equation (4.9), with the result 


B=pJ /2ak. 


Fig. 4.2 A general path of integration around a long straight wire. 


From symmetry, the vector B must be perpendicular to the wire and to the line OP. For a small 
movement d/ along the path the contribution to the line integral is 


B-dl= BR do =! ao (4.10) 
TW 


This expression is independent of the radius and so, for any closed path encircling the wire, 
f B-dl=y,/ (4.11) 


This result is known as the magnetic circuit law. It may be stated in words as follows: The line 
integral of the vector B around a closed path is equal to [lo times the total current enclosed. This form 
of the law is derived from Equation (4.11) by using the principle of superposition. An equation in the 
form of (4.11) is written for the magnetic flux density due to each conductor and then the set of 
equations is summed. The result remains true when the current is distributed as a current density 
rather than being confined in a set of wires. 


The general form of the law can be written 


(4.12) 


f B-di=y,|[ J-da 
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where the normal component of the current density J is integrated over an open surface S which spans 
the closed curve as illustrated by the figure on the right. This equation, like Gauss’ theorem in 
electrostatics, can be applied usefully to problems only when it is possible to deduce the distribution 
of the magnetic flux density from the symmetry of the problem. We shall see in Chapter 5 that it has a 
very important use in the solution of problems in which the greater part of the magnetic flux is carried 
by an iron circuit. 


4.4 Magnetic scalar potential 


Although the theory of magnetism has been developed here without using the idea of magnetic poles, 
it is useful to introduce the magnetic analogue of electric potential. This concept is a convenient 
mathematical device which does not have any physical significance. The magnetic scalar potential 
difference between the points P and Q, by analogy with Equation (1.13), is 


T 


1 oe 
o-t,=-—| B-dl (4.13) 
Mo 


where the constant /é9 has been introduced so that Equation (4.13) has the same form as the equations 
given in older texts which assume the existence of magnetic poles. 
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The magnetic scalar potential is of limited usefulness 

because it can be given a unique value only ina 

region of space in which there are no electric p 

currents. This point is illustrated by the figure on the 

right. It can be shown from Equation (4.10) that the Z 
difference in the scalar potential between P and Q is 

independent of the path of integration chosen 54) ‘ 
provided that it does not encircle any electric current. 

If, however, a path such as 2 were chosen for the 

integration, the result would be to subtract / from the 

value of the potential difference calculated. 


What is worse still is that, by making several loops around the current, it is possible to make the 
potential difference take a whole series of different values, because the value of the integral changes 
by J for each complete loop made around the current. Thus the magnetic scalar potential can be given 
a unique value only in a region of space containing no currents. Despite this limitation there are still 
times when the idea is useful. 


The differential form of Equation (4.13) is 


B=-y,grad ? =—-“,VT (4.14) 


Moreover, since it has been assumed that free magnetic poles do not exist, the magnetic analogue of 
Gauss’ theorem is 


{} B-dA=0 (4.15) 


or, in words, the flux of B out of a closed surface is zero. The unit of magnetic flux is the weber (Wb). 
The differential form of Equation (4.15) is 


divB=V-B=0 (4.16) 


Combining Equations (4.14) and (4.16), we find that must obey Laplace’s equation in free space 
just as the electric potential does. That is 


div(grad # )=V°# =0 (4.17) 


The commonest use of the magnetic scalar potential is for calculating the magnetic flux density in the 
space between iron pole pieces. We shall return to this subject in Chapter 5. 
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4.5 Forces on current-carrying conductors 


We have already seen that an element of a current-carrying conductor can be regarded as equivalent to 
a moving charge. Substitution from Equation (4.5) into Equation (4.4) gives the force on a current 
element as 


F =/(d1,B) (4.18) 


The total force on a circuit carrying an electric current can be found by integrating Equation (4.18) 
over all the current elements in the circuit. This is of central importance in electric power engineering 
but rarely needed in electronics. You should consult other texts for details of its application. 


4.6 Summary 


The law of force between two electric charges moving with constant velocities differs from the 
electrostatic law of force by the presence of an additional term. This additional component of force is 
the magnetic force. Magnetism is thus shown to be a particular aspect of electricity rather than a 
separate phenomenon. The concept of magnetic flux density is a useful tool for calculating the 
magnetic forces on charges and currents. There are two ways of calculating the magnetic flux density 
produced by electric currents: the Biot-Savart law; and the magnetic circuit law. 


The concept of scalar potential, which is useful in electrostatics, is found to be less useful in 
magnetostatics. It can be given a unique definition only in a region of space which is free from 
current-carrying conductors. The main application of the magnetic scalar potential is in calculating the 
magnetic flux density distribution in air gaps in magnetic circuits. The question of the magnetic force 
on a current-carrying conductor has been touched on only briefly because of its relative lack of 
importance in electronics. 
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5. The magnetic effects of iron 


Objectives 


To show that the polarization of iron in a magnetic field can be represented by an equivalent 
surface current. 

To introduce the magnetic field vector H as a means of calculating magnetic fields in the 
presence of iron without needing to calculate the equivalent surface current distribution. 

To introduce a revised form of the magnetic circuit law which is valid in the presence of 
magnetic materials. 

To derive the boundary conditions at the interface between two materials with different 
magnetic properties. 

To show that iron can be regarded as a conductor of magnetic flux and to consider the 
application of the idea to magnetic screening. 

To discuss the calculation of fields produced by current-carrying conductors in the presence 
of iron by the method of images. 

To apply the analogy between electric current and magnetic flux to magnetic circuits, 
introducing the concepts of magnetomotive force and reluctance. 

To consider approximate ways of dealing with fringing and leakage fluxes in magnetic 
circuits. 

To discuss the description of the properties of a magnetic material by its hysteresis loop, 
including the concepts of remanence and coercive force and the distinction between hard and 
soft magnetic materials. 
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e To consider the solution of magnetic circuit problems when non-linear effects are important. 

e To consider the solution of magnetic circuit problems involving permanent magnets, 
including questions of the stabilization of the magnets and Evershed’s criterion for optimizing 
the design of a circuit. 


5.1 Introduction 


All engineering students must, at some time, have played with magnets and been fascinated by their 
ability to pick up, not just single iron objects, but whole strings of pins or paperclips. The strange 
property of iron and one or two other materials, which we call ferromagnetism, is employed in many 
ways in electronics. 

e The ability of ferromagnetic materials to concentrate and conduct magnetic flux is used in 
inductors and recording heads. 

e Permanent magnets find application in loudspeakers. 

e Special high-resistivity materials known as ferrites are used to make aerials, inductors for 
radios and microwave components such as circulators and isolators. 

e Magnetic media are used for recording purposes in the form of audio and video tapes and 
computer disks though this is less important than it used to be. 

e The most important uses of ferromagnetic materials, if sheer weight of material were to be the 
criterion, are in electrical power engineering, in transformers, generators and motors. These 
lie outside the scope of this book, but they are important examples of the application of the 
fundamental ideas in electromagnetism discussed in this chapter. 


5.2 Ferromagnetic materials 


The overwhelming majority of materials which exist in the world have magnetic properties which 
differ so little from those of a vacuum that, for nearly all engineering purposes, they are 
indistinguishable from it. The exceptions to this rule are the elements iron, cobalt and nickel, together 
with certain materials which include them. Their special properties are illustrated by Fig. 5.1. 


Figure 5.1(a) shows an air-cored solenoid carrying a current /. The flux density in the vicinity of the 
solenoid, somewhat as shown, could be mapped with the aid of a flux meter. It is worth emphasizing 
that there would be no measurable change in the field if the solenoid were filled with a core made of 
brass, aluminium, wood, plastic, or any other non-ferromagnetic materials. If, on the other hand, an 
iron core were inserted the result would be a considerable increase in the flux density, as shown in 
Fig. 5.1(b). 
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( p i 
ee (a) 


{a) 
(b) 
Fig. 5.1 The magnetic fields of (a) an air-cored solenoid and (b) a similar solenoid with an iron core, 
showing diagrammatically the great increase in the strength of the field in the second case. 


This property of iron arises because the iron atoms are themselves tiny 


permanent sources of magnetic flux like current loops. If each of these ‘loops’ is i {7 

represented by an arrow to show its polarity, then a qualitative explanation of a 4 
ferromagnetism can be given. In an un-magnetized piece of iron the atomic Unmagnetized 
magnets are arranged, head to tail, in closed loops. All the flux then circulates =: = 
within the iron and none of it outside. The effect of applying an external flux —. 


density to the iron is to cause the atomic magnets to line up with the external flux Megnetized 


so that their flux is added to it. These two situations are illustrated 
diagrammatically in the figures on the right. The actual processes of 
magnetization are more complex than this, but there is no room in this book to 
enter into them in more detail. The macroscopic phenomena of ferromagnetism 
are discussed later in this chapter. 


From Fig. 5.1(b) we observe that the field outside the solenoid could just as well be produced by the 
addition of a second winding, with a suitable current and distribution of turns, to the air-cored 
solenoid of Fig. 5.1(a). This suggests that, just as the polarization of dielectric materials can be 
represented by a surface-charge distribution, so the polarization of ferromagnetic materials can be 
represented by a surface-current distribution. This is consistent with the view adopted in this book that 
magnetic phenomena are best regarded as resulting from the motion of electric charges. 
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NV turns 
of wire 


Fig. 5.2 An iron ring of mean radius R with N turns of wire wound uniformly upon it. 


To illustrate the effect of the magnetization of iron we consider an iron ring with N turns of wire 
wound uniformly upon it as shown in Fig. 5.2. Let the current in the winding be J, and let the 
magnetization of the iron be represented by an additional current /’ in the winding. Then, by the 
magnetic circuit law (Equation (4.11)) the flux density in the ring is obtained from 


22RB = N(I+1') (5.1) 
or 
N 
ere Oey 5.2 
sap toll +1’) (5.2) 
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Now experiments with electromagnets show that the magnetic flux density produced increases with 


the current, though not necessarily in a linear fashion, so let us assume that J’ = vy, J where is 
g y Kon Lai 


the magnetic susceptibility. 


Then 
N 
B=— pll+ ee) 
rR’ ( Ln) ( ) 
or, if fs, = 14743 
N N 
| 7 gece 54 


where the quantities 4 and “,, are known as the permeability and the relative permeability of the 
iron, respectively. The analogy with the case of dielectric permittivity is only superficial because it is 
hardly ever valid to regard the magnetic permeability as a constant, The other distinguishing feature is 
that the relative permeability of iron is much greater than the relative permittivity of common 
materials, being of the order of 10°. 


In problems involving dielectric materials we found it useful to make use of the vector D to avoid 
having to calculate the distribution of polarization charges. The same approach can be used in 
magnetic problems. We define a new vector H by the equation 


HB (5.5) 
mW 


so that, for the simple case shown in Fig. 5.2, we have 


NI 


H =—— 5.6 
27R me) 


The vector H is known as the magnetic field. It is measured in A m’ in SI units, as can be seen from 
Equation (5.6). Since the strength of the field is proportional to the product of the current in the 
winding and the number of turns, it is common to speak of a field strength of so many ‘ampere-turns 
per metre’. Equation (5.6) also shows that the vector H depends only upon the distribution of electric 
currents and not upon the arrangement of any iron present. 


Making use of the vector H, we can deduce a form of the magnetic circuit law (Equation (4.12)) 
which is valid when magnetic materials are present: 
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f H-di={[J-da 


To find the distribution of magnetic flux density in problems involving 
iron we first calculate H from the known distribution of currents, using 
Equation (5.7), and then calculate B in each region from Equation (5.5). 


Magnetic problems cannot usually be solved by the application of potential theory and the use of 
Laplace’s equation because of the non-linearity of magnetic materials. For the same reason it is not 
possible to make use of the principle of superposition unless linearizing assumptions can be made. 
This contrasts with the situation in electrostatics where the distribution of charges is frequently 
unknown but the potentials are given on the conducting boundaries. 


5.3 Boundary conditions 


Whenever a magnetic problem involves two or more regions of space with different magnetic 
properties it is necessary to use the relationships between the magnetic vectors on the two sides of 
each boundary. 


The figure on the right shows the boundary between two materials having permeabilities “2, and /,. 


Consider a Gaussian surface in the form of a ‘pillbox’ as shown. If the cross-sectional area of the box 
in the plane of the surface is dA, then the flux of B out of the box is 


(B,,—B,,)dA=0 (5.8) 


from Equation (4.15), since the contribution of the flux through the sides of the box becomes 
negligible as the height of the box tends to zero. The consequence of Equation (5.8) may be stated in 
words: The normal component of B is continuous across a boundary. 


The boundary condition on H is found by applying the magnetic circuit law to 
a small closed path encircling part of the boundary as shown in the figure on 
the right. If the length of the path parallel to the boundary is dl, then the line 
integral of H around it becomes 


(H,, —H,, )dl = current enclosed (5.9) 


It should be noted that the current referred to is an actual electric current. It does not include the 
hypothetical surface currents which were introduced as a convenient way of describing magnetic 
materials. In many cases the surface current will be zero and then the boundary condition is that the 
tangential component of H is continuous across the boundary. 
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5.4 Flux conduction and magnetic screening 


Equation (5.9) has one very important practical corollary. Consider an iron bar surrounded by air. 
Assuming that there is no surface current flowing along the bar, we know that the tangential 
components of H are equal inside and outside the bar. But, since B = “H, the component of the 


magnetic flux density parallel to the bar is much greater inside it than outside it. 


If such a bar is placed in a magnetic field in which the total flux is 
fixed, with the axis of the bar aligned with the field, then the i eee 
requirement that B;,.,>> B,i, means that most of the flux will pass 

through the iron as shown in the figure on the right. The flux density 

in the air is then much less than it was before the iron was introduced. i 


__ een OO 


a ee 


The iron therefore acts as quite a good conductor of magnetic flux. 


The analogy between flux conduction by iron and electric conduction is strengthened still further if we 
consider what happens to the direction of a flux line entering the iron from the air, as shown in Fig. 5.3. The 
boundary conditions require that B,; = Bj, and that B,; = ,B,. Now 1, is typically in the range 10° - 10°, so the 
flux line must turn through an angle even greater than that indicated in Fig. 5.3. The angle of incidence is 
greatest when the flux line within the iron lies nearly parallel to the surface. A simple calculation shows that if 
a = 89° and = 1000, then a, = 3.3°. For many purposes it is adequate to make the approximation that iron 
is a perfect magnetic conductor, so that all the flux lines enter it at right angles to the surface. 
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Fig. 5.3 The change in direction of a magnetic flux line when passing from air into iron. 


The ability of iron to conduct flux is put to work in a number of ways. It allows magnetic circuits to 
be constructed which conduct the flux wherever it is required - an idea which we shall consider in 
detail later on. It also allows flux to be excluded from regions where it is not required. Figure 5.4 
shows what happens when a hollow iron cylinder is placed in a magnetic field. The flux lines prefer to 
pass through the iron and the space within the cylinder is almost entirely free from magnetic flux. 
This property can be used to shield sensitive equipment from magnetic fields, and special alloys such 
as mumetal are used for this purpose. 


lron cylinder 


Fig. 5.4 An iron cylinder acts as a conductor of magnetic flux and screens the volume inside it from 
the external flux. 


The effectiveness of a screen can be expressed in decibels as 


(5.10) 


B at a point in the absence of the screen 
B at the same point with the screen 


Sy = 20 oe 


A mumetal cylinder 100 mm in diameter and | mm thick has a magnetic screening effectiveness of 
60dB provided that the external magnetic field is not strong enough to saturate it. It should be noted 
that it is only low-frequency fields which are screened out by the mechanism of flux conduction. A 
different mechanism operates at radio frequencies. 


The fact that iron acts as a good conductor of flux means that it is often possible to calculate the 
distribution of flux around iron pole pieces by regarding them as equipotential surfaces and solving 
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Laplace’s equation (4.17). This method can only give relative values of the flux density because the 
potentials of the pole faces are not known. The absolute values can be found by using the fact that the 
flux crossing the air gap is equal to that in the iron, apart from any leakage flux, and applying the 
magnetic circuit law to the whole flux path. The application of the magnetic circuit law to problems of 
this type is discussed in the next section. 


The fact that an iron surface can be regarded as a magnetic equipotential surface allows the method of 
images to be used for the solution of problems in magnetostatics. 


5.5 Magnetic circuits 


The analogy between the conduction of electricity by conductors and the conduction of magnetic flux 
by iron leads us to the useful concept of a magnetic circuit. Figure 5.5 shows a simple magnetic 
circuit formed by a square iron core with a narrow air gap in it at A. A coil of N turns of wire is wound 
on the core and carries current /. To make the problem easy to handle we make some simplifying 
assumptions. Let us suppose that the magnetic field strength is the same everywhere within the iron, 
having the value H;, and that the magnetic field in the air gap is H,. 


Air gap 


Fig. 5.5 A simple magnetic circuit made up of an iron core with an air gap in it. The flux is supplied by 
a winding of N turns of wire. 


Now consider a closed path around the circuit which follows the centre line (shown dotted). Applying 
the magnetic circuit law to this path gives 


4LH,+gH, = NI (5.11) 


Strictly speaking, the length of the iron path is (4Z - g) but, provided that g is small, the difference is 
negligible given the other assumptions which have been made. 


The boundary condition (5.8) requires that 


(5.12) 
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where ® is the total flux circulating, and A is the cross-sectional area of the iron at right angles to the 
direction of the magnetic field. We have therefore assumed that all the flux due to the coil is contained 
within the bounds of the iron core and its projection across the air gap. We also know from Equation 
(5.5) that 


and BL =MH, (5.13) 


Substitution for H; and H, in Equation (5.11) in terms of ® and the various constants gives 


nr=t [es £g (5.14) 
ALM My 


which can be written 
oe = 00 (5.15) 


where {| (= N/) is known as the magneto-motive force and as the reluctance. 
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Equation (5.15) is analogous to Ohm’s law for electric circuits, but the analogy must not be pressed 
too far. Unlike an electric current, the magnetic flux in the circuit dissipates no energy and there are 
no circulating magnetic charges. Moreover, it cannot be emphasized too often that magnetic materials 
do not normally behave in a linear fashion, so the reluctance of a circuit can only be regarded as 
approximately constant at best, and that only under a limited range of conditions. We will return to 
give this point further consideration after discussing the behaviour of real magnetic materials. First, 
though, it is useful to mention two other matters: fringing and leakage of magnetic flux. 


5.6 Fringing and leakage 


In the discussion of magnetic circuits above it was assumed that the flux lines 
in the air gap passed straight across, as shown in the upper figure on the right. 


A little thought shows that this is an over simplification and that the flux 
lines must actually spread out into a fringing field, as shown in the lower 
figure in the margin. The truth of this observation can be demonstrated by 
laying a sheet of paper over the air gap and sprinkling iron filings on it. 
Clearly the reluctance of the magnetic circuit is affected by the fringing of 
the field around the gap so we need an estimate of how big the effect is. 


A detailed solution of the problem would require the solution of Laplace’s equation in the air gap 
(assuming the iron surfaces to be equipotentials). This is unnecessarily complicated for most 
purposes. A useful rule of thumb is that the effective cross-sectional area of the air gap is larger than 
that of the pole face by a strip of width g/2 around it, as shown in Fig. 5.6. Equation (5.12) then 
becomes 


O=Bw =B(w+g) (5.16) 


L 


and the reluctance of the circuit is 


o-| ca -| (5.17) 


Fig. 5.6 The fringing of the magnetic field in an air gap can be allowed for, approximately, by adding 
an extra strip around the area of the pole face. 
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Fringing fields are put to use in magnetic recording heads. Figure 5.7 shows a typical arrangement. 
Most of the flux in the magnetic circuit passes directly across the air gap, but some of the fringing 
field passes through the magnetic coating on the tape and magnetizes it. To achieve linear 
performance the signal supplied to the coil is superimposed on an a.c. bias current whose frequency is 
at least five times the maximum signal frequency. 


The tape is replayed by passing it across either the same head or a similar read head where the varying 
magnetization of the tape produces a varying flux in the magnetic circuit. We shall see in Chapter 6 
that varying the flux passing through a coil causes an electromotive force to be induced in it. The 
e.m.f. can be processed to recover the original signal waveform. 


Coil 
Recording head 
Air gap 
Magnetic coating 


Tape —— 


Tape backing 
Fig. 5.7 The arrangement of a magnetic recording head. 
A second difficulty in calculating the flux in the air gap is illustrated by Fig. 5.8. Some of the flux 


simply takes a short cut such as A-A’ and never reaches the gap at all. This flux is known as leakage 
flux. It exists because, although iron conducts flux better than air, the air is a poor magnetic insulator. 


Fig. 5.8 An electromagnet showing a typical flux leakage path A—-A’. 


If we wish to pursue the analogy with electric circuits it is necessary to think of the material 
surrounding the conductor as being a poor conductor (such as wet earth) rather than an insulator. 
Fortunately, in most cases, a very crude estimate of the leakage flux is adequate to provide a 
satisfactory value for the flux in the air gap. 
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5.7 Hysteresis 


In the preceding sections the permeability of iron has been taken to be a constant. While this is a 
useful first approximation, which can be justified in a limited number of cases, it is not generally 
possible to regard permeability as a constant. The relationship between the magneto-motive force 
applied to a specimen of iron and the magnetic flux density produced is a far from simple one. The 
magnetization of the iron depends not just on the present value of the magneto-motive force, but also 
on the previous history of the specimen. In this book we are restricting our attention to macroscopic 
phenomena, so no attempt is made to offer an explanation of their causes in terms of the atomic theory 
of matter. 


One way of investigating the relationship between the magneto-motive force applied to a specimen of 
iron and the flux density produced would be to use a magnetic circuit like the one shown in Fig. 5.5. 
Provided that the air gap is very small, it is possible to neglect the effects of fringing. The flux density 
in the air gap is then the same as that in the iron and its strength could be measured by placing a Hall 
probe in the gap. The total magneto-motive force applied to the circuit is N/, so the magneto-motive 
force applied to the iron is 


. 2 B 
= N-@ = n=. (5.18) 
Ho 
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Now |; _ is the product of H; and the length of the part of the circuit lying within the iron, so it 
would be possible to deduce values of H; from measurements of the current in the coil and the flux 
density in the air gap. The results of such measurements are customarily displayed by plotting B; 
against H;. There are better ways of obtaining the information than the experiment described, but we 
are concerned here with the properties of the iron rather than the details of how they are measured. 


Fig. 5.9 A typical hysteresis loop showing the initial magnetisation path (ab), saturation (at b), 
remanence (at c), the coercive force (at d), and an inner loop (shown dotted). 


Figure 5.9 shows a typical B—H curve. If the specimen of iron is initially un-magnetized, then its state 

can be represented by the point a at the origin of the graph. Now suppose that the current in the coil is 
gradually increased. It is found that the B-H curve follows the path ab, rising slowly at first, then more 
rapidly, and, finally, levelling off. This curve is known as the initial magnetization curve of the iron. 
When the curve levels off at high values of H; the material is said to be saturated. 


If the current in the coil is reduced to zero, the flux density does not retrace the curve ab. Instead, it 
follows a curve such as bc, so that the flux density in the circuit does not fall to zero when the current 
is zero. The value of B; at c is referred to as the remanence of the iron. 


In order to reduce the flux in the circuit to zero it is necessary to reverse the current in the coil. The 
value of H; needed to make B; zero (the point d@) is known as the coercive force. Increasing the current 
in the coil beyond the level needed to reach d eventually produces saturation of the iron at e with the 
direction of the flux opposite to that at b. Finally, reducing the current to zero and then increasing it 
with reversed polarity produces the curve efgb. 


If the current in the coil is repeatedly taken through the same cycle, the hysteresis loop bcdefgb is 
traversed repeatedly in a stable manner, provided only that the maximum value of the current in the 
coil is the same for both polarities. In many applications the flux in the iron is produced by a periodic 
(though not necessarily sinusoidal) current, and the iron then behaves in the manner described. If the 
maximum current is increased still further, there is no change to the loop because the magnetization of 
the iron cannot be increased beyond its value at saturation. If, on the other hand, a periodic current of 
smaller amplitude is used, it is found that the behaviour of the iron is described by a smaller loop such 
as the one shown dotted in Fig. 5.9. 
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The properties of magnetic materials can be varied by making alloys with different proportions. They 
are also strongly influenced by the ways in which the finished material is prepared, especially any 
heat treatments used. A very great variety of materials now exists. Their properties can be looked up 
in reference books. From a practical point of view the main division among magnetic materials is 
between ‘hard’ and ‘soft’ materials. The differences between these two classes are illustrated in Fig. 
5.10. 


Soft 


Hard 


Fig. 5.10 Typical hysteresis loops for hard and soft magnetic materials. 


Soft magnetic materials are characterized by having narrow hysteresis loops, low remanence and 
small coercive forces. They are therefore easily magnetized and demagnetized and are used as 
conductors of flux in magnetic circuits and magnetic screens. Because these materials have narrow 
hysteresis loops their behaviour can be approximated by their initial magnetization curves. Moreover, 
if the loop is also nearly straight for fields below saturation, like the one shown in Fig. 5.9, it is 
possible to make the approximation that the permeability is constant. 


Hard magnetic materials have broad hysteresis loops, high remanence and high coercive forces. 
Such materials are difficult to demagnetize. They are used for making permanent magnets and 
magnetic recording materials. 


A further important feature of the properties of magnetic materials is illustrated by Fig. 5.11. Suppose 
that the material has been magnetized to saturation (b) and then demagnetized to the point p by the 
application of a reverse field. If the magnetizing field is then reduced to zero, the working point of the 
material moves to g along a minor loop. Minor loops are usually narrow and it is possible to 
approximate the loop between p and q by the straight line pq, which is usually nearly parallel to be. 
Provided that the demagnetizing field applied does not take the working point past p back onto the 
main loop, the material will operate in a stable manner along the minor loop pq. This fact is of 
considerable importance in the application of permanent magnets. 
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Fig. 5.11 Part of a B—H curve showing a minor loop p-q. 


5.8 Solution of problems in which p cannot be regarded as constant 


Consider again the magnetic circuit shown in Fig. 5.5. Neglecting leakage and fringing fluxes we may 
assume that B;= B,. The application of the magnetic circuit law to the circuit gives 


B, 
NI =4LH, + 9H, =4LH, +24 

Ho 
or 


B, = “°(NI-4LH,) (5.19) 
g 
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A second condition is that B; and H;must be related to each other by the hysteresis loop of the 
material: 


8, = f(H,) (520) 


Equations (5.19) and (5.20) are a pair of simultaneous equations in B;and H;. However, because (5.20) 
is not the equation of a straight line, they are non-linear simultaneous equations. 


They can be solved graphically, or by using an analytical approximation to the equation of the 
hysteresis loop. The graphical solution is shown in Fig. 5.12. The working point Q is at the 
intersection of the straight line represented by Equation (5.19) with the hysteresis loop. The intercept 
of this line on the horizontal axis is the total magneto-motive force in the circuit, while its slope is 
determined by the relative sizes of the iron path and the air gap. The method is similar to the process 
of finding the working point of a transistor by drawing a load line across the characteristic curves. 


NI 
4L 


Fig. 5.12 A magnetic circuit problem can be solved by drawing a load line across the B-H plot for the 
iron. The working point of the circuit is at Q. 


5.9 Permanent magnets 


We have already noted that some magnetic materials can be magnetized so that they produce a 
substantial magnetic flux even when the magnetizing field is removed. These hard magnetic materials 
are used for making permanent magnets. Permanent magnets are used in loudspeakers and in a 
variety of microwave devices such as ferrite isolators, magnetron oscillators and travelling-wave tube 
amplifiers. 


Figure 5.13 shows a typical B—H plot for a permanent magnet material. This curve, which is the part 
of the hysteresis loop lying in the second quadrant, is known as the demagnetization curve of the 
material. 
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oO 


Fig. 5.13 The working point for a permanent magnet circuit lies in the second quadrant of the B—H 
plot at a point such as Q. A permanent magnet can be stabilized against demagnetization by external 
fields by operating it at a point Q’on a minor loop such as R-S. 


Suppose that the circuit shown in Fig. 5.5 is made of this material and magnetized to saturation by 
passing a current through the coil with the air gap bridged by a piece of soft iron. When the current in 
the coil is reduced to zero the working point will be at P in Fig. 5.13. Neglecting fringing and leakage 
fluxes we have, when the soft iron is removed from the gap, 


4L 
B, = (bot (5.21) 
& 


from Equation (5.19), setting /= 0. This load line passes through the origin and has negative slope. It 
is plotted in Fig. 5.13 as the line OQ. The effect of opening the gap is therefore to demagnetize the 
magnet to some extent. For this reason it is common to speak of the demagnetizing field of the air 


gap. 


If the piece of soft iron is reinserted in the air gap, the operating point of the magnet moves along a 
minor loop to the point P’. In this state the magnet is largely immune to the effects of external fields. 
For this reason it is usual to store permanent magnets with their air gaps bridged by a piece of soft 
iron known as a ‘keeper’. 


Suppose that someone has forgotten to put a keeper on the magnet and that it is exposed to a magnetic 
field which tends to demagnetize it. The working point might then move to R. When the field is 
removed the working point lies at 0’ on the minor loop RS. This illustration shows that it is necessary 
to treat permanent magnets carefully if their properties are not to be affected. For example, it is 
unwise to try to force a pair of magnets together with their fields opposing because they are then 
trying to demagnetize each other. 
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It is not a good thing that the working point of the magnet should be sensitive to external influences in 
this way. For this reason permanent magnets are often stabilized by deliberately demagnetizing them 
beyond their working points on the main hysteresis loop. The operation of magnets which have been 
stabilized is stable along a minor loop (or recoil line) such as RS unless the external demagnetizing 
field is very strong. Smaller external fields produce a temporary shift in the working point, but it 
returns to QO’ when the perturbing field is removed. 


5.10 Using permanent magnets efficiently 


In a circuit involving a permanent magnet, the cost of the magnet is usually a considerable part of the 
cost of the whole circuit. It is, therefore, desirable to use the magnet material as efficiently as possible. 
A tule of thumb for this purpose can be derived by considering the circuit shown in Fig. 5.5. Let the 
lengths of the magnet and the air gap be /; and /, and let their cross-sectional areas be A; and A,. The 
magnetic circuit law gives 


LH,+lA,=0 (5.22) 
and the conservation of flux gives 


A,B, = A,B, (5.23) 
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The volume of the magnet is 


V, =14, =-1,4, 2 | 


iB. 


L L 


Now /,A, is the volume of the air gap which may be taken as given, and H, and B, are the field 
strength and flux density in the air gap. So, the minimum magnet volume to provide a given flux 
density in a given air gap is achieved by making the product H;B; as great as possible. This is known 
as Evershed’s criterion. The quantity HB; is known as the energy product of the magnet for reasons 
which will become apparent in a later chapter. 


Figure 5.14 shows the curves of B against H and BH against B as they are usually plotted. The 
optimum working point of the material is then that which gives the value of B corresponding to the 
maximum value of BH as shown at Qpp. 


Fig. 5.14 Typical curves of H and BH versus B for a permanent magnet material showing the choice 
of working point at the maximum value of BH to make the most efficient use of the magnet. 


5.11 Summary 


In this chapter we have considered the properties of ferromagnetic materials. We have seen that their 
principal characteristic is that they acquire a very strong magnetization when placed in a magnetic field. As 
with dielectric materials it is possible to represent the properties of the materials by introducing a new vector, 
in this case H, which includes the effects of the polarization of the materials. Using the vector H a new form 
of the magnetic circuit law was derived, which is valid even when magnetic materials are present. 


Consideration of the boundary conditions at the surface of a ferromagnetic material showed that such 
materials can usually be regarded as good conductors of magnetic flux. This property allows the path 
of the flux to be controlled so that it is directed to the place where it is needed, as in a recording head, 
or directed away from sensitive areas in magnetic screening. 


In a limited range of cases it is possible to regard ferromagnetic materials as being linear, but, in 
general, it is necessary to take account of their non-linear behaviour as described by their hysteresis 
loops. 
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6. Electromagnetic induction 


Objectives 

e To show that an electromotive force is induced in a conductor moving through a magnetic 
field, and that a current flows if the circuit is completed. 

e To show that an electromotive force is induced in a loop of wire moving through a non- 
uniform magnetic field and that, by implication, an e.m.f. is induced in a loop of wire when 
the magnetic flux linked to it changes. 

e To generalize the first two objectives in the form of Faraday’s law of electromagnetic 
induction. 

e To demonstrate the links between the circuit concepts of self and mutual inductance and 
electromagnetic field theory. 

e To discuss the causes of electromagnetic interference and ways of reducing It. 

e To introduce methods, including energy methods, for calculating self and mutual inductances. 

e To introduce the analogy between L, C and R in two-dimensional problems. 

e To discuss the concept of the storage of energy in a magnetic field, and to show the 
equivalence between calculations of the energy stored in an inductor from field and circuit 
points of view. 

e To discuss the special case of energy storage in a magnetic field in iron and to introduce the 
idea of hysteresis loss. 

e To introduce the induction of eddy currents in a conductor by a changing magnetic field. 


6.1 Introduction 


In this chapter we turn our attention to phenomena involving conductors and magnetic fields in 
relative motion, and magnetic fields which are changing with time. These effects, first investigated by 
Faraday in 1831, complete the link between magnetism and electricity. We have already seen how 
magnetic fields are produced by electric currents. We shall now show that electromotive forces and 
electric currents can be produced by electromagnetic induction. These ideas are fundamental to the 
generation of electricity by electromagnetic machines and to the use of transformers in the distribution 
of electric power. As far as electronic engineers are concerned their importance lies in the fact that 
they provide an explanation of self and mutual inductance, and the means to calculate these circuit 
parameters. 


6.2 The current induced in a conductor moving through a steady 
magnetic field 


Consider a piece of straight wire moving with constant velocity v through a 

steady magnetic field B which is at right angles to the wire and to its direction 

of motion, as shown in the figure on the right. The conduction electrons in the e 
wire are carried through the field with the wire so, from Equation (4.4), each 
electron experiences a magnetic force gvB acting upon it. The cross product 
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symbols have been omitted because the directions of the wire, the velocity of the wire and the field 
have been chosen to be mutually perpendicular. As a result of this force the electrons are displaced so 
that the ends of the wire acquire polarization charges, as shown in the figure. 

The polarization of the wire produces an electric field within it which opposes the motion of the 
electrons. In equilibrium the electric force on each charge must be equal and opposite to the magnetic 
force, so that 


E=vB (6.1) 


The magnitude of the potential difference between the ends of the wire is therefore 
V= | Edl=vBl (6.2) 


Now suppose that the ends of the wire are connected to parallel rails by sliding contacts and that the 
circuit is completed by a load resistor, as shown in Fig. 6.1. The potential difference between the ends 
of the rod acts as an electromotive force in the circuit and a steady current flows such that 


IR=vBl (6.3) 
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The electric power dissipated in the resistor is 

?R=vBII (6.4) 
The mechanical power input is the product of the force needed to move the wire and its velocity. 
Since the current in the circuit is /, the force needed to move the wire is B//, from equation (4.18), so 
the mechanical power input is 


Fv=vBlI (6.5) 


This is exactly equal to the electric power, demonstrating that this device converts mechanical energy 


LE 


Fig. 6.1 A current is induced in a closed circuit in a magnetic field if the shape of the circuit is 
changed. 


into electrical energy. 


Examination of Fig. 6.1 shows that the force exerted on the current-carrying slider by the field is in 
the opposite direction to v. This internal force must be opposed by an equal and opposite external 
force if the wire is in steady motion. The external force is in the same direction as v, so it does work 
on the slider, providing a power input to the device. This simple arrangement forms the starting point 
for the discussion of the generation of electricity by electromagnetic machines. Important as the topic 
is, it is outside the scope of this book and will not be pursued any further here. 


Recalling the definition of magnetic flux which was used in earlier chapters, we can see that, with the 
direction of velocity shown, the flux linked to the circuit of Fig. 6.1 is changing at a rate given by 


a® __wp (6.6) 
dt 


but, from Equation (6.2), this is just the electromotive force developed in the circuit. The direction of 
the electromotive force is such that when a current flows in the circuit the magnetic forces oppose the 
external mechanical force which is causing the flux to change. These ideas can be combined by 
writing 


a (6.7) 
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6.3 The current induced in a loop of wire moving through a non- 
uniform magnetic field 


In the experiment described in the previous section it was shown that a current could be induced in a 
circuit by altering the shape of the circuit and so changing the flux linked to the circuit. It is 
interesting to see whether a current could also be induced by keeping the shape of the circuit fixed and 
varying the strength of the magnetic field. This could be achieved by taking a square loop of wire, as 
shown in Fig. 6.2(a), and moving it with a steady velocity v through a magnetic field perpendicular to 
the plane of the loop whose variation with position in the direction of motion is given by the graph of 
Fig. 6.2(b). 


(a) (b) 
Fig. 6.2 A current is induced in a closed loop of wire when it is moved through a non-uniform 
magnetic field at right angles to the plane of the loop. 


At the instant when the wire AD is at x relative to the origin, the e.m.f. generated in it is 
x 
%,,, =vIB, =) (6.8) 


where the order of the subscripts shows that the polarity of the e.m-f. is such as to drive a current from 
A to D through the rest of the circuit. In the same way the e.m.f. generated in BC is 


x+l 


=, =vlB, (2+) (6.9) 


The remaining sides of the loop are moving parallel to their own directions so no potential differences 
are generated between their ends. The net e.m-f. in the circuit is therefore 


e= via =) =v/B, (6.10) 


directed so as to tend to drive a current in a counter-clockwise direction. Now the flux linked to the 
circuit is 
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o=ra 1?) (6.11) 
so that 

d® dx 

—=/B,—=vIB 6.12 
dt 0h VID ( ) 


So, in this case also, the e.m.f. is equal to the rate of change of flux linkage. The direction of the flux 
produced by the induced current is opposed to the original field, so that it tends to prevent the flux 
linked to the loop from increasing as the loop is moved. Since it is the relative motion of the loop and 
the field which produces the e.m.f. it is reasonable to suppose that the same result would have been 
produced if the loop had been held fixed and the source of the field moved with velocity v in the 


negative x-direction. 


6.4 Faraday’s law of electromagnetic induction 


In the previous sections it has been shown that, at least in the two cases considered, an electromotive 
force is induced in a circuit when the magnetic flux linked to the circuit is changed. This phenomenon 
was observed experimentally by Faraday, who generalized his findings in the law which bears his 
name. 
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Faraday’s law of electromagnetic induction states that, if the flux linked to a circuit is changed in any 
way, then an electromotive force is induced in the circuit whose magnitude is proportional to the rate 
of change of the flux linkage to the circuit. 


The term flux linkage is used here to include circuits such as coils in which the wire encircles the 
path of the circuit more than once. In these cases the flux linked to the circuit is the sum of the flux 
linkage to the individual turns. It is often possible to make the approximation that the flux linkage to 
each turn of a coil is the same, so that the total flux linkage is the product of the flux linked to one 
turn and the number of turns. Denoting flux linkage by 4, this relationship is expressed by the 
equation 


A=N® (6.13) 


We have also seen that the direction of the e.m.f. induced in a circuit by a changing magnetic flux 
linkage is always such that it tries to oppose the change of flux which causes it. This statement is 
known as Lenz’s law. 


Faraday’s law can be expressed in mathematical terms by writing 
o= —— (6.14) 


For a simple circuit such as that shown in Fig. 6.3, the flux linkage is given by 


A=[B-dA (6.15) 
S 


where the integration is carried out over an open surface spanning the circuit and the elementary 
vectors dA are taken to point upwards as shown. 


Be at 


induced 


rr = 


Fig. 6.3 The direction of the current induced in a loop of wire by a magnetic flux which is increasing 
with time. 
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The e.m.f. in the circuit can be calculated in a similar manner by making use of Equation (3.17). The 
direction of the line integral around the circuit is chosen to be in a right-hand corkscrew sense with 
respect to the direction of the vector dA. This direction is shown in Fig. 6.3 by the direction of the 
elementary vector dl. Thus 


of E-dl (6.16) 


Now, if B is in the direction shown and it is increasing in magnitude with time, Lenz’s law requires 
the direction of the induced current to be opposite to the direction of the line integral as shown in Fig. 
63: 


Substituting for ! and 4 in Equation (6.14) produces a general mathematical form of Faraday’s law: 
fe-d=-“{/B-4a (6.17) 
dt 


Although this expression may seem abstract and complicated it should be remembered that it is no 
more than a restatement of Faraday’s law and Lenz’s law in the symbolic language of mathematics. 


When a circuit has more than one turn, the integral on the right-hand side of Equation (6.17) is taken 
over all the turns. In the problems which are of interest in electronics the circuit is normally fixed and 
the magnetic field changing with time. Under these circumstances the differentiation can be brought 
inside the integral to give 


éB 
p E-dl=—|{—-da (6.18) 


6.5 Inductance 


We are now ina position to consider the third of the trio of passive circuit parameters, resistance, 
capacitance and inductance. Figure 6.4 shows two loops of wire arranged so that when a current is 
passed through loop | some of the flux produced passes through loop 2 and vice versa. 


ee 8B 
i= > 
1 
Fig. 6.4 Two loops of wire are inductively linked when some or all of the magnetic flux produced by 
one of them passes through the other. 
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Suppose that the flux linked to loop 2 is 42; when the current in loop | is /;. If the loops are in a region 
of constant permeability, then 42; is proportional to /,. The constant of proportionality is 


Gb. (6.19) 


where ZL, is the mutual inductance coupling the loops. From Equation (6.19) it is evident that 
inductance has the dimensions weber per ampere, and this unit is called the Henry (H). The ratio of 
the flux linked to loop | to the current in the same loop is the self-inductance of that loop, given by 


L,,=—t (6.20) 


Two similar equations can be derived by assuming that the source current is in loop 2. 
These equations can then be brought together by using the principle of superposition to give 


A, =Ay +Aj, HLL, + Ly!) 
A, Sh the Blof tlel, (6.21) 
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As this is a passive linear system the matrix of the coefficients of inductance must be symmetrical; 
that is 


L,, =L,=M (6.22) 


where M is the mutual inductance of the loops. Using this notation the double subscripts can be 
dropped and the self inductances represented by the symbols L; and L>. When Faraday’s law is applied 
to equations (6.21), we obtain 


weap thy ee 
dt dt 
dl, _ dl 
® =--M—!_-f, — 6.23 
, yo a (6.23) 


where the positive directions of the currents and the e.m.f.s are the same. These are the familiar circuit 
equations for inductively coupled circuits. It is important to remember that these equations hold only 
when the system is linear, that is, when the fluxes are directly proportional to the currents producing 
them. Since the behaviour of iron can be highly non-linear, as we saw in Chapter 5, it cannot be taken 
for granted that iron-cored inductors behave as linear circuit elements. 


The self and mutual inductances in a circuit are related to each other, as will now be demonstrated. 
Consider again the pair of loops shown in Fig. 6.4, but suppose that loop | is a coil of N; turns and 


coil 2 has N2 turns. When a current J; flows in | a flux ®, is generated. For the sake of simplicity it 
will be assumed that the whole of this flux is linked to every turn in coil 1, so that 


A, =N,®, (6.24) 
In general, not all of the flux will pass through coil 2 so that the flux linked to that coil can be written 
A, =k,N,®, (6.25) 


where k is the proportion of the flux generated by coil 1 which is linked to coil 2. It follows 
thatk, <1. 


Assuming linearity, we can write 


M A, _&k,N, 
= = (6.26) 
L, A, N, 
If a current is passed through coil 2 we obtain, by a similar argument 
M A, _ &KN, 
= = (6.27) 
L, A, N, 
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where k; is the fraction of the flux produced by coil 2 which passes through coil 1. Multiplying 
Equations (6.26) and (6.27) together gives 


Mahi, =r Ld, (6.28) 
where the coupling coefficient k =./k,k, <1. 


6.6 Electromagnetic interference 


In Chapter 2 we saw that there can be unwanted capacitive coupling between two circuits. Similar 
problems can arise through inductive coupling. Figure 6.5 shows a simple example of the coupling of 
two circuits through their mutual inductance M. If an alternating current flows in circuit I, then it 
produces a magnetic flux, some of which may be linked to circuit 2. 


Fig. 6.5 When two circuits are coupled by stray mutual inductance the signal in one circuit can 
interfere with that in the other. For example, if circuit 1 is the a.c. mains a spurious 50 Hz signal could 
be added to the signal V2 in circuit 2. 


To simplify matters we will assume that /, is much greater than /, so that, approximately, 


1 
I, =——— JV exp|jo t 
1 R,, +R, 1 pj it) 


and the e.m.f. induced in circuit 2 is 


‘0 M 
V.=-Mo1=-L°"" 7 exp(jo t) (6.29) 
dt Ry + Ry 


This voltage is added to the signal voltage in circuit 2. The interference is worst when V; is large and 
the loop impedance of circuit | is small. The commonest example of this is when circuit | represents 
the a.c. mains and the interference appears as 50 Hz ‘mains hum’. In this case it is not possible to 
reduce the interference by altering V; or (Rs; + Rz1) and effort has to be concentrated on reducing M. 
To do this we have to reduce the proportion of the flux produced by circuit | which is linked into 
circuit 2. 
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Three strategies are possible: 

1. Reduce the area of circuit 2 
If the signal source V> is a transducer some distance from the amplifier, then this can be 
achieved by using a pair of wires twisted together (a ‘twisted pair’) to connect them together. 

2. Rotate circuit 2 so that its plane is parallel to the flux produced by circuit I 
Rotating one circuit relative to the other can reduce the flux linkage to zero. If the two circuits 
are to be enclosed in the same case in fixed positions this can be a useful approach. 

3. Puta screen of high permeability magnetic material in a position which screens circuit 2 
from the flux of circuit 1. 
The special alloy known as mumetal (1, ~ 10°) is used for this purpose. This method works 
only at low frequencies. 


A particularly troublesome type of electromagnetic interference is caused by earth loops. Figure 6.6 
shows a typical situation in which this problem occurs. Two electronic instruments, A and B, are 
connected to the main supply by three-core cables whose earth conductors are connected to the 
instrument frames. The instruments are also connected to each other by a coaxial cable whose sheath 
is connected to the frames of the instruments to screen the signal wire from capacitively coupled 
interference. This arrangement produces a closed earth loop as shown by the cross-hatching in Fig. 
6.6. 
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~ 5 
Mains 


— R 


Fig. 6.6 An earth loop, shown shaded, can be formed by the earth conductors of the electricity supply 
to two electronic instruments and the sheath of the coaxial cable connecting them together. 


The earth loop has low resistance and large currents can be induced in it by the flux of the mains or of 
some other source of interference. The induced current (/;) flows through the resistance of the sheath 
of the coaxial cable (R-) as shown in Fig. 6.7. This produces a spurious potential difference between P 
and QO which is added to the signal voltage Vy. 


Fig. 6.7 The potential difference across the resistance Rc of the sheath of the coaxial cable in Fig. 6.6 
produced by the circulating earth loop current /; is in series with the signal source Va. 


Two solutions to this problem are possible: 
1. Break the earth loop so that no current can flow in it. 

This is not as simple as it seems. Disconnecting the earth lead from the plug at R or S is potentially 
very dangerous because lethal voltages could appear on the instrument frames under fault conditions. 
This solution is only possible if a permanent earth connection between A and B can be ensured, for 
example by bolting them to the same rack. Ideally there should only be one path to earth from any 
point on the equipment. The alternative of disconnecting the cable sheath at P means that the signal 
earth follows the path PSRO. The mutual inductance between this circuit and the mains is high. 


2. Increase the resistance of the earth loop. 
This can be achieved by putting resistors between the signal earth and the frame of each 
instrument as shown in Fig. 6.8. If the added resistors (Rp) are much larger than Rc then the 
circulating current in the earth loop is reduced because the e.m.f. is fixed. The unwanted 
voltage appearing across Rc is reduced by a factor which is approximately Rc/2Rp. It is 
important that Rp is small enough so that no part of the system which can be touched reaches 
a lethal voltage when the maximum fault current is flowing. If Rc = 0.1 Q and Rp = 50 Q, 
then the interference voltage is reduced by a factor of 1000. If, also, the mains fuse of each 
instrument is rated at 250 mA, then the voltage at P or Q cannot exceed 12.5 V with respect to 
earth. 
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Fig. 6.8 The resistance of the earth loop can be increased by connecting resistors Rp between the 
signal earth and the supply earth of electronic instruments. 


6.7 Calculation of inductance 


Before proceeding to a discussion of the ways of calculating inductance it is useful to review the 
situations in which such a calculation might be necessary. They are: 
e To find the self and mutual inductances of circuit components such as inductors and 
transformers. 
e To find the inductance per unit length of two-wire transmission lines. 
e To estimate the stray mutual inductance between parts of a circuit such as parallel tracks on a 
printed circuit board. 


Inductance is calculated either by assuming a distribution of currents and computing the flux linkage, 
or by assuming a flux distribution and then finding the ampere turns needed to produce it. The first 
approach is best for many purposes, but we shall see later in the chapter that the two methods lead to 
upper and lower bounds for the inductance when energy methods are used. The calculation of 
inductance in simple cases may be illustrated by two examples. 


Figure 6.9 shows a simple transformer in which the primary and secondary windings having NV, turns 


and N, turns respectively are wound upon an iron core whose permeability is w. The mean path length 
of the magnetic circuit is L and its cross-sectional area is A. 


Fig. 6.9 The arrangement of a simple transformer. 
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Let a current /; flow in the primary winding, then, by the magnetic circuit law the magnetic field in 


the core is 
H = Ni, 
L 


and the flux in the core is 


® =BA=pHA= I 


1 


UNA 
L 


The fluxes linked to the primary and secondary windings are 


2 
A NA 
i ey and FG ee acl cas 
L 
so that the self and mutual inductances are 
2 
A ,A 
|e ie ace and M= 2 _ HNN, 
f, L I, L 
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In the same way by assuming that a current /> flows in the secondary winding we obtain 


UN>A A yah HNN 
i ie i 


A 
L,=—= 6.34 
2-7 (6.34) 


We have shown, incidentally, that Equation (6.22) is valid for this case and that MZ = L,L, as would 


be expected from Equation (6.28) for the perfect coupling between the windings which has been 
assumed here. 


In real transformers there is always some leakage of flux, so that the inductances calculated in this 
way are only approximate. To achieve greater accuracy it would be necessary to compute the detailed 
distribution of flux in and around the core. The effects of leakage can be minimized by using 
multilayer windings. 


At frequencies above 50 Hz the capacitance between the windings can be important because it limits 
the band of frequencies over which the transformer will work. The complete transformer can be 
represented by an equivalent circuit and methods exist for calculating the parameters of the circuit for 
any particular design of transformer. The core can be regarded as linear to a first approximation if it is 
made of a soft magnetic material and operated well below saturation. Inductors (‘chokes’) only have a 
single winding. They are sometimes made with an air gap in the magnetic circuit, so that the 
reluctance of the circuit is dominated by the reluctance of the gap and is, therefore, linear. 


Figure 6.10 shows unit length of a coaxial cable whose inner and outer conductors have radii a and b 
respectively. We wish to find an expression for the series inductance per unit length of the line. 


Fig. 6.10 A coaxial cable 


If the ends of the cable are short-circuited then the two conductors form a closed circuit as shown. If 
the magnetic circuit law is applied to a circular path which is concentric with the cable then net 
current is only enclosed when the path lies in the space between the conductors. Thus there is a 
magnetic field only in the space between the conductors and, from considerations of symmetry, it 
must be directed in the tangential direction as shown. 
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Let the current in the conductors be / then using the magnetic circuit law we find that the magnetic 
field strength on a path between the conductors at radius r is 


ie = (6.35) 
2nr 


The conductors form a closed circuit having only a single turn, so the flux linked to it is found by 
calculating the flux crossing a unit length of a radial plane, as shown in Fig. 6.10. Thus, assuming that 
the material in the space between the conductors is non-magnetic 


a2nr a 


b I 
A= i‘ PO alt (2) (6.36) 
The inductance per unit length is the flux linkage per unit current 


i oo ®) (6.37) 
I 2a a 


The significance of this result will become apparent in Chapter 7 when the use of coaxial cables as 
transmission lines is discussed. The estimation of the inductance per unit length of transmission lines 
of other shapes can be achieved by energy methods. This topic is treated later in this chapter. 


6.8 Energy storage in the magnetic field 


In Chapter 2 we saw that the work done in setting up an electric field could be thought of as being 
stored either in a lumped manner in the capacitance of the system ('ZCV°) or distributed throughout 
space with an energy density '2D.E. In just the same way the work done in setting up a magnetic field 
can be regarded as being stored either in the circuit inductance or distributed throughout the magnetic 
field. 


It is well known from elementary circuit theory that the energy stored in an inductor is given by 


ae ere 
W =| Li di =5H (6.38) 
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An alternative expression for W is obtained by substituting A// for L 


W =5Al (6.39) 


To derive an expression for the energy density in a magnetic field 
we consider a current loop and a typical flux tube linked to it, as 
shown in the figure on the right. The stored energy associated 
with the tube, obtained from Equation (6.39), is 


AW = 5 Ab! (6.40) 


Now, if the cross-sectional area of the flux tube at a particular point is 6A and the flux density at the 
same point is B, then 


A® =B-OA (6.41) 
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Applying the magnetic circuit law to the flux tube gives 
[=$H-dl~> Hol (6.42) 


if the flux tube is considered to be made up of volume elements of length 6/ as shown in the figure 
above. Substituting in Equation (6.40) for A® and J produces 


AW = +(BoA)y. (H61)= >) (BH6A61) (6.43) 


where the vector dot products have been dropped because the vectors concerned are all parallel to 
each other from the definition of the flux tube. (B 6A) is a constant of the flux tube, so it is permissible 
to bring it inside the summation. Finally, noting that the product (dA6/) is just the volume, dv, of the 
element defined by them, the energy stored in the flux tube becomes 


AW = D5 BH dv (6.44) 


so that the energy density in the magnetic field is evidently ‘2BH. As in the case of the electric field 
considered in Chapter 2, a more rigorous derivation shows that the energy density in the field is given 
in general by 


~ SB-H (6.45) 


so that the total energy stored may be written 


W =—[[JB-Hav (6.46) 


106 


Electromagnetism for Electronic Engineers 6. Electromagnetic induction 


6.9 Calculation of inductance by energy methods 


We have seen that energy methods provide a useful way of obtaining quite good approximations to 
capacitance and resistance in cases which cannot be solved by analytically. These methods can also be 
used to estimate inductances. The argument runs exactly parallel to the cases of capacitance and 
resistance so it will not be given in detail. Starting from Equation (6.38) we discover that 


(6.47) 


where W'is the estimate of the stored energy produced by assuming that the current distribution is 
given. The corresponding lower bound can be obtained by making use of the definition of self- 
inductance to give another expression for the stored energy 


a, NX’ 


f= 
21, 


(6.48) 


If an estimate W” of the stored energy is obtained by assuming that the flux distribution is given, then 


6.49 
2w" es 


6.10 The LCRZ analogy 
From Equation (2.10) it is easy to show that the shunt capacitance per unit length of a coaxial cable is 


2 
CaF (6.50) 


In(b/a) 


A comparison between Equations (6.37) and (6.50) shows a marked resemblance. This suggests that, 
if we can calculate the capacitance for a particular two-dimensional arrangement of electrodes, we 
should be able to deduce an inductance for the same geometry from it. The justification for this is that 
the field patterns in each case must satisfy Laplace’s equation (1.27). Since this also applies to current 
flowing in a conductor, we can add resistance to the discussion. We will assume that the medium 
between the electrodes has the same properties everywhere. 


It is easy to show that the resistance per unit length for radial current flow between concentric 
electrodes (Fig. 6.11(a)) separated by a uniform conducting material of resistivity p is given by 


R= In(b/a) (6.51) 
20 
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6. Electromagnetic induction 


Fig. 6.11 (a) Coaxial and (b) parallel strip geometries 


Similarly, for the parallel strip geometry shown in Fig. 6.11(b) we obtain the following parameters per 
unit length if fringing fields are neglected: 


C=ew/d (6.52) 
L=yud/w (6.53) 
R= pd/w (6.54) 
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We deduce that for all such two-dimensional arrangements of electrodes the parameters per unit 
length are given by 


C=eT (6.55) 
Lepr (6.56) 
R=p/F (6.57) 


where the parameter I is determined solely by the geometry of the electrodes. For coaxial geometry 


2 

f= (6.58) 
In(b/a) 
and for parallel strips 
(6.59) 

T=w/d 
By multiplying Equations (6.55) and (6.56) together we obtain 
LC=eu (6.60) 


which is a constant whose value depends only on the properties of the material surrounding the 
electrodes. 


Dividing Equation (6.56) by Equation (6.55) and taking the square root we get 


E-f4()-2 (6.61) 
C é\T T 


where Z has the dimensions of resistance and depends only upon the properties of the material. 


Equations (6.55) to (6.57) can be used to find any two of the three parameters C, L and R for a particular 
geometry once one of them is known. Note carefully that the magnetic flux lines in each case are at right 
angles to the electric field lines. The significance of this will become apparent in Chapter 7. 


These ideas can be extended further by invoking the principle of duality discussed in section 3.6 in 
which the roles of the flux lines and equipotentials are exchanged. It is easy to see from Fig. 6.11(b) 
that the parameters of the dual system are given by 


C'=e/T (6.62) 
L' =u (6.63) 
R=. (6.64) 
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6.11 Energy storage in iron 


So far our discussion of energy in magnetic fields has been restricted to those cases where the 
permeability is constant. To investigate storage of energy when the permeability is not constant, we 
consider an iron ring of mean circumference / and cross-sectional area A which has N turns of wire 
wound upon it, as shown in Fig. 6.12. 


Fig. 6.12 The arrangement of a toroidal inductor 


In time dt let the current in the winding increase from / to (/ + d/) and the flux in the iron increase 
from @ to (® + d®). The induced e.m.f. in the winding is NV d®/ dt so the work done in the time 


interval is 


aW =NIPadt=NId0 (6.65) 
But 

@®=BA and HI=NI 

so 

dW = AlH dB (6.66) 


The change in the stored energy when the system is taken from flux density B, to By is given by 

By 
W=Al| HH dB (6.67) 
But A/ is just the volume of the iron, so the change of energy density is 


By 
we I, H dB (6.68) 
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Once again more rigorous argument shows that Equation (6.68) is correct for all cases if the scalar 
product H dB is replaced by a vector dot product H- dB . The work done in magnetizing an initially 
unmagnetized specimen of iron is illustrated by Fig. 6.13. The integral in Equation (6.68) is 
represented by the shaded area between the initial magnetization curve and the vertical axis. 

8 


Fig. 6.13 The work done in magnetizing a piece of iron is represented by the area between the initial 
magnetization curve and the vertical axis. 
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6.12 Hysteresis loss 


A case of particular interest is the work done in taking a piece of iron once around its hysteresis loop. 
In any iron-cored inductor or transformer the core is taken through this cycle for every cycle of the 
a.c. current in the windings. Figure 6.14 shows a typical hysteresis loop. The change in stored energy 
per unit volume in going from | to 2 is represented by the shaded area in Fig. 6.15(a). Let us denote 
this energy by wi. 


Fig. 6.14 The work done in taking a piece of iron around its hysteresis loop is proportional to the area 
of the loop. 


Fig. 6.15 (a) The increase in the stored energy of a piece of iron when it is taken around its hysteresis 
loop from 1 to 2; (b) the energy recovered when it is taken from 2 to 3. 


The shaded area in Fig. 6.15(b) represents the change in stored energy per unit volume when the iron 
is taken from 3 to 2, a change which is mathematically possible but physically impossible. This 
energy 1S W32. The real physical change is obtained by exchanging the limits of integration in Equation 
(6.68) so that w23 = -w32. The change in energy density resulting from going from | to 3 via 2 is 


Wi3 = Win — W39 


which is represented by the difference between the areas shown in Figs 6.15(a) and 6.15(b), and by the 
shaded area in Fig. 6.14. Since the hysteresis loop is symmetrical the work done in encircling the loop 
once is equal to the area of the loop multiplied by the volume of the iron. It is evident that this work is 
not zero since all real ferromagnetic materials have hysteresis loops with non-zero areas. 
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This phenomenon means that electrical energy is changed into heat in the iron core at a rate 
proportional to the frequency of the signal applied to the device. The conversion of energy in this way 
is known as hysteresis loss. The loss is minimized by using soft magnetic materials for transformer 
cores, as an examination of Fig. 5.10 will show. 


6.13 Eddy currents 


When any conducting circuit is placed in a changing magnetic field currents are induced in the circuit, 
as we saw earlier. This remains true when the circuit is a solid piece of metal. In that case currents, 
known as eddy currents, circulate within the metal, causing power loss by ohmic heating. They are 
put to use in the industrial eddy-current heating (ECH) process, which is used for brazing metal 
components together. In other cases, particularly in transformer cores, they give rise to unwanted 
losses. 


The full theory of the generation of eddy currents is beyond the scope of this book, but it is possible to 
give an approximate treatment of the special case in which the conductor is in the form of a thin strip 
with the magnetic field parallel to its length. 


Fig. 6.16 The element of volume used in estimating the eddy-current loss in a thin conducting strip. 


Figure 6.16 shows a section of a strip of conducting material whose width (/) is much greater than its 
thickness (2f). The strip is placed in a uniform magnetic field directed along the strip with flux density 


By cos(at). The eddy currents in the strip flow in closed loops like the one shown. From Faraday’s 


law we know that the e.m.f. induced in the loop is given approximately by 
dB : 
ea ery = 2) x@B, sin(at) (6.69) 
since the area of the loop is approximately 2/x. 
The resistance of the loop is approximately 
2pl 
dx 


per unit length of the strip, where p is the resistivity of the material. 


R= (6.70) 


The mean power loss in the loop is 
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a aT (Ba 
fa 


so that the mean power loss per unit length of the strip is 


2p2 ; 3,2 p2 
jee) ae eo (6.71) 
a al 3p 


We can also express this as mean power loss per unit surface area 


3,2 p2 
W, = ra By (6.72) 
3p 


In deriving this expression we have implicitly assumed that the eddy currents are not strong enough to 
have a significant effect upon the strength of the magnetic field. 
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Equation (6.72) shows that the eddy-current losses depend very strongly on the thickness of the strip. 
This is the reason why transformer cores are made of thin strips or laminations of steel which are 
insulated from each other by a coating of lacquer. They are normally made of special steel which has 
a high resistivity because, as can be seen from Equation (6.72), this also helps to reduce the losses. 
The mean power loss also increases rapidly with increasing frequency. At radio-frequencies the losses 
in laminated iron cores are unacceptably high and magnetic oxides of iron known as ferrites are used 
instead because they have much higher resistivities. 


6.14 Real electronic components 


By this stage it should be apparent to the reader that the behaviour of real electronic components is 
much more complicated than that of the lumped resistors, capacitors and inductors used in circuit 
theory. Whenever current flows through a component a magnetic field is generated, so resistors and 
capacitors must have some self-inductance. The resistance in the winding of an inductor or 
transformer means that there must be some voltage drop across the device even under d.c. conditions. 
The existence of a potential difference between the different parts of the winding implies that there is 
some capacitance present. Although these parasitic effects can usually be neglected at low 
frequencies, they must be included in the circuit representations of the components at high 
frequencies. 


6.15 Summary 


In this chapter Faraday’s law of electromagnetic induction has been introduced. It has been shown 
that an electromotive force can be induced in a circuit which has magnetic flux linked to it either by 
changing the area of the circuit or by changing the strength of the magnetic field. From this starting 
point the idea of inductance was introduced, so establishing another link between electric circuit 
parameters and electromagnetic field theory. 


Stray inductances between parts of electronic circuits can produce unwanted coupling. The causes of 
electromagnetic interference were considered, and some cures suggested. The calculation of 
inductance was discussed using both direct calculation and energy methods. The latter are based on 
the idea of the stored energy associated with a magnetic field being distributed throughout the field. 
The analogy between L, C and R in two-dimensional problems was introduced as a way of deducing 
these parameters from each other. 


The special case of energy storage in iron was discussed, leading to an expression for the work done 
in taking a sample of iron once around its hysteresis loop and to the idea of hysteresis loss. Finally, 
the consideration of the effects of a changing magnetic field on a solid conductor led to a discussion 
of eddy currents and of the losses associated with them. 
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7. Transmission lines 


Objectives 

e To introduce electromagnetic waves through a discussion of transmission lines as distributed 
circuits. 

e To discuss the solutions to the wave equation for sinusoidal waves and pulses, and to 
introduce the idea of phase velocity. 

e To explain the concept of characteristic impedance and its relationship to the power flow in 
the line. 

e To explain the use of complex notation to describe waves. 

e To show that a wave is reflected at the end of a line unless it is terminated by a load equal to 
its characteristic impedance. 

e To introduce the terms ‘voltage reflection coefficient’ and ‘voltage standing wave ratio’. 

e To show how impedances are transformed by a transmission line, and to introduce the idea of 
the quarter-wave transformer as a matching device. 

e To consider the field description of a transmission line and to show how the phase velocity 
and characteristic impedance can be calculated. 

e To derive expressions for the electric and magnetic fields in a coaxial cable and to consider 
how the power flow in the line may be calculated from them. 


7.1 Introduction 


Two-wire transmission lines appeared early in the history of the practical application of electricity 
with the advent of the electric telegraph. They are still an important means of transmitting both 
electrical power and information, although for the latter purpose they are now supplemented by radio 
and microwave links and optical fibres. The commonest forms of transmission line are parallel-wire 
and twisted-pair lines, used for electricity distribution and telephone connections, and coaxial cables, 
used for television aerial downleads and for interconnecting electronic instruments. 


Throughout this book I have tried to show the links between the field and the circuit approaches to the 
description of electromagnetic phenomena. These provide alternative ways of dealing with problems. 
Electronic engineers are generally happiest when they are able to use circuit methods, with the 
properties of the components represented by equivalent circuits. Field theory provides an alternative 
to experimental measurements as a means of determining the equivalent circuit parameters. In this 
chapter, however, the field and circuit approaches draw even closer together. We shall see that the 
transmission of signals on transmission lines can be described in terms of propagation of 
electromagnetic waves. We shall also see that, when the dimensions of a circuit become comparable 
with the wavelength of these waves, elementary circuit theory breaks down. 
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7.2 The circuit theory of transmission lines 


It was shown in Chapters 2 and 6 that a two-wire transmission line such as a coaxial cable has shunt 
capacitance and series inductance uniformly distributed along its length. At low frequencies it is 
usually possible to ignore these impedances and treat circuits as sets of lumped components connected 
together by wires whose length has no effect upon the operation of the circuit. We must now 
investigate the circumstances in which this approximation is not valid and the consequences for the 
design of electronic circuits. 


/ Ye Ldx Yeldx /+d/ 


V Cdx VidV 


Fig. 7.1 The 'tee' network representing a short length dx of a transmission line having series 
inductance L per unit length and shunt capacitance C per unit length. 


Figure 7.1 shows the circuit representation of an element of transmission line of length dx. If the line 
has inductance L per unit length and capacitance C per unit length, then the inductance and 
capacitance of the element must be L dx and C dx. 
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The inductance has been divided into two so that the network is a symmetrical tee section. For the 
moment we assume that the whole line stretches to infinity in either direction, being made up of 
identical tee sections joined together in a cascade. If the voltage across the capacitor is V’, then 


Pap pee 0 ea ea 
2 dt 2 dt 


When the second-order term ad is neglected we obtain 
ov — as (7.1) 
Ox Ot 


The current through the capacitor is 


dl = Gg 
dt 


To first order V'=V , so that 


ol OV 
Sc cdl (7.2) 
ox Ot 

The pair of simultaneous differential equations (7.1) and (7.2) defines the relationship between the 
voltage and the current in the limit when the length (dx) of the element shrinks to zero and the line 
becomes completely uniform. 


The current can be eliminated from these equations by differentiating Equation (7.1) with respect to x, 
and Equation (7.2) with respect to ¢, since 


Or al 

Oxot otox 

giving 

ov ov 

ae =LG ae (7.3) 


Similarly, differentiating Equation (7.1) with respect to ¢ and Equation (7.2) with respect to x, we get 


Orr O71 
—— SiC 74 
Ox? Ot? ee 
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Equations (7.3) and (7.4) are examples of the one-dimensional form of the wave equation. For 
sinusoidal waves their general solution has the form 


V =V, cos(@t +kx)+V, sin(ot + kx) (7.5) 
and 
I =1,cos(at +kx)+ I, sin(ot + kx) (7.6) 


where V;, V2, J; and J, are constants whose values are determined by the boundary conditions of the 
problem. Substituting these expressions into the wave equations shows that they are acceptable 
solutions provided that 


k? =@’° LC CD) 


The physical significance of these solutions can be discussed, without loss of generality, by 
considering the behaviour of the function 


V =V, cos(@t — kx) (7.8) 


The form of this function is shown in Fig. 7.2. The solid curve shows the variation of voltage with 
position along the line when t = 0. At a later time ¢ = ¢; the wave has moved to the position shown by 
the broken line. This can be verified by noting that the first zero of the function is given by 


(ot —kx) = 


none " 


This equation shows that the zero crossing, and hence the whole wave, is moving in the positive x- 


N/a 


direction with a constant velocity given by 
a) 

y=— (710) 
k 


This velocity is known as the phase velocity of the wave. 


119 


Electromagnetism for Electronic Engineers 7. Transmission lines 


Fig. 7.2 The variation of the voltage along a transmission line, as given by Equation (7.8) for t= 0 and 
a later time t = t; showing how the wave moves steadily in the positive x-direction. 


From Equation (7.7) we also have 


oe (7.11) 


Equation (6.60) shows that v, depends only on the properties of the material between the conductors 
and not on their geometry. 
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The propagation constant, or wavenumber, & is related to the wavelength A of the wave by 
k =— (7.12) 


as can be seen from Fig. 7.2. If kis negative the phase velocity is also negative and the wave is 
travelling in the negative x-direction. This corresponds to the positive signs in Equation (7.5). The 
alternative solution, involving the sine function in place of the cosine, gives a wave which is shifted in 
phase by 90°. 


By making a suitable choice of the constants V; and V2 in Equation (7.5) we can make the amplitude 
and phase of the wave what we will. This is made clearer by a consideration of the alternative form of 


(7.5): 


V =V,cos(@t+kx+ g) (7.13) 


where the amplitude Vo and the phase @ are arbitrary constants whose values are to be determined by 


the boundary conditions. The equivalence of the two forms of the wave given in Equations (7.5) and 
(7.13) can be demonstrated by expanding the cosine function in Equation (7.13). 


7.3 Representation of waves using complex numbers 

In the theory of a.c. circuits ‘j-notation’ is an indispensable tool. Before proceeding to a discussion of 
the properties of transmission lines it is necessary to extend the notation to problems where the 
voltages and currents vary sinusoidally in distance as well as in time. To do this consider the voltage 


V' = jV, sin(@t —kx)— jV, cos(at — kx) (7.14) 


When this voltage is added to the general solution for the wave in the positive x-direction, obtained 
from Equation (7.5) by taking the negative signs, the result is 


V =V,|cos(wt —kx)+ jsin(wt — kx)|— /V, [cos(at — kx)+ jsin(ot —kx)| 
or 


V =Vexp|j(at — kx)| (7.15) 


where the amplitude of the wave V is, in general, complex. It turns out that Equation (7.15) provides 
a more convenient form of representation for a wave than Equation (7.5), but it must be remembered 
that it is only the real part of Equation (7.15) which has physical significance. 
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Hence, whenever the form in Equation (7.15) is used to represent a wave in a problem it is implicit 
that, when an expression has been obtained as the solution to the problem, the real part of it is to be 
taken as having physical significance. It is easy to verify by substitution that Equation (7.15) is a 
solution of Equation (7.3). The two arbitrary constants in Equation (7.5) are present as the real and 
imaginary parts of V. 


7.4 Characteristic impedance 


So far we have discussed the propagation of waves on the line solely in terms of the voltage. The 
corresponding solution for the current can be obtained by substituting 


V=V, exp| j(ot — kx)| and I=], exp| j(at — kx)] 
into Equation (7.1), with the result 
— fkV, =—joLl, 


so that, making use of (7.7) 


L= 2 pa (7.16) 
oL Ly 
where 
oL iL 
Z. =-—=_JO— TAT 
aa C (7.17) 


This is a constant which has the dimensions of resistance. It is known as the characteristic 
impedance of the line. Equation (6.61) shows that Z) depends upon the geometry of the line and the 
properties of the material between the conductors. For the loss-less line considered here Zp is a real 
quantity, that is, a pure resistance, because the voltage and current are in phase with each other. It can 
be shown that if, on the other hand, the line has loss, then the voltage and current are no longer in 
phase with each other and, consequently, Zp) has an inductive or capacitive component. The physical 
significance of the characteristic impedance becomes clearer when we consider the reflection of 
waves from the termination of the line. 
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7.5 Reflection of waves at the end of a line 


Figure 7.3 shows a transmission line of characteristic impedance Zp which is terminated by an 
impedance Z at x = 0. 


—————+  Vrexp j(ot - kx) 


es 


V,.exp j(at + kx) 


x=0 
Fig. 7.3 A line of characteristic impedance Z) terminated by an impedance Z with incident and 
reflected waves having amplitudes V; and V,. 


The wave incident on the termination is 
V =V,exp| j(ot —kx)| (7.18) 


We must assume, until it has been shown otherwise, that some of the incident wave will be reflected 
as the wave 


V =V, exp|j(ot + kx)| (7.19) 
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The amplitudes V; and V, are complex. The corresponding currents are 


l= fi Jest j(ot —kx)| (7.20) 
Zo 
and 
i= { ie Jesh j(ot + kx) (7.21) 
Zo 


as can be shown by substituting the expressions for the incident and reflected voltage waves in 
Equation (7.16). We shall see later that the negative sign in Equation (7.21) indicates that the power in 
the reflected wave is travelling in the negative x-direction. 


At x = 0 the voltage is 


V= (V, aoe Jexp(jat) (7.22) 


L Tr 


and the current is 


[= e l Jes jat) (7.23) 


0 
by adding together the voltages and currents for the incident and reflected waves. 
But the ratio of the voltage to the current at x = 0 must equal the termination impedance Z. Therefore 


V,+V, 


Z= Zz, (7.24) 


which can be rearranged to give 


Ve DAD 
V. Z+Z, 


(7.25) 


where I is known as the voltage reflection coefficient. Equation (7.25) shows that the amplitude of 
the reflected wave is zero when Z = Zp. When this condition is satisfied the termination is said to be 
matched to the line. 


When a line is terminated by a load equal to its characteristic impedance the instantaneous power 
absorbed in the termination is 
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11 | cos*(at) (7.26) 


The mean power is obtained by averaging Equation (7.26) over one complete cycle, giving 


se 727 
TA (7.27) 


Since there is no reflected wave when Z = Z) it follows that W is the mean power flow in the incident 
wave. This is the condition for maximum transfer of power from the incident wave to the load and it 
is a generalization of the Maximum Power Transfer Theorem of circuit theory. 


If the load resistor were replaced by a semi-infinite section of transmission line having the same 
characteristic impedance there would still be no reflected wave in the part of the line to the left of x = 
0. Thus the characteristic impedance is the input impedance of the semi-infinite line. Since we could 
just as well have replaced the load resistor by a section of line of finite length terminated by Zp 
without producing a reflected wave it follows that the input impedance of a transmission line 
terminated by its characteristic impedance is also Zo. 


When the line is not terminated by its characteristic impedance the voltage at a general point on it is 
V =V, exp j(ot —kx)]+ TV, exp| (at + kx) (7.28) 
The amplitude of the voltage which would be measured at this point is 


|V|=|V,| exp(— jx)||1+ TD exp(2 jkx)| =|V;||1+ Dexp(2 jkx)| (7.29) 


since | exp(- jkx)| =1. 


Equation (7.29) shows that the amplitude of the voltage varies along the line between the maximum 
and minimum values 


| 


max 


y= 


max 


V,|(1+ ||) (7.30) 
V;|(1-|C|) (7.31) 


The ratio of these two is easily measured with a sliding probe and a suitable detector. It is known as 
the voltage standing wave ratio (VSWR), and its value is given by 
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HT 
Ar 


(7 32) 


The probe, therefore, detects a signal whose amplitude varies periodically with position along the line 
as a standing wave. The wavelength of the standing wave observed is given by 2kx = 27 from 
Equation (7.29), that is 


x= 


4 (7.33) 
5 ; 


where A is the wavelength of the travelling waves. 


7.6 Pulses on transmission lines 
So far we have assumed that the signals propagating on the transmission lines are sinusoidal. It is, of 
course, possible to consider non-sinusoidal signals by using Fourier synthesis to construct them from 


sine waves of different frequencies. But pulses, which are important for digital data transmission, can 
most easily be studied directly. It is easy to show by substitution that any functions having the form 


V =V, f(x—v,t) and T=1,f(x-v,t) (7.34) 


are solutions of Equations (7.1) to (7.4). 
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The function f can have any form, but for the study of 
pulses it is taken to be the step function shown in the figure 
on the right. This function is defined to be 0 for 


(x—v,t)> 0 and 1 for (x—v,t)< 0. Substituting the 


voltage and current given in Equation (7.34) into Equation 
(7. 1) gives 


Vif' =v Lif’ 

where f” is the derivative of fso that, making use of Equations (7.11) and (7.17) 

Vy =Zoly (7.35) 
as before. 


Figure 7.4(a) shows a voltage source with internal impedance Zs; connected to a semi-infinite 
transmission line of characteristic impedance Zp. If the line is uncharged until the switch is closed at 

t =0, then a wave having the form shown in Fig. 7.4(b) will propagate down the line with velocity v,. 
As v, is independent of frequency the step propagates without any change in its shape. A line for 
which this is true is said to be non-dispersive. Conversely, if v, depends on frequency, the line is 
dispersive, the different Fourier components of the step travel at different speeds, and the shape of the 
step changes as it propagates. 


Zs 
Vs Zo 
(a) 
V 
—____—__» 
VsZo / (Zo Gs Zs) 
(b) x 


Fig. 7.4 A semi-finite transmission line of characteristic impedance Z, is connected to a source of 
impedance Zs by a switch at t= 0. The result is a step voltage wave which travels down the 
transmission line as shown. 
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The voltage and current at the start of the line must satisfy Equation (7.35), so that 


Zo 


ef 7 36 
ve veo 


0 
by the potential-divider rule. 


7.7 Reflection of pulses at the end of a line 


When a transmission line is terminated by a resistance R, part of an incident pulse is usually reflected. 
The voltage and current at any point on the line are given by 


V =V,+V, and (take 
Z) 

To satisfy the boundary conditions at the termination 
pattie 7 

lf ¥,=¥, 
so that 
a a 
= (7.37) 
Vv, BAZ, 


From Equation (7.37) it is clear that there are three possible conditions: 


1. R> Zp: V,is positive so that the reflected wave is added to the incident wave. 
2. R<Z: V,is negative and the reflected wave is subtracted from the incident wave. 
3. R=Z): V,1is zero and the incident wave is completely absorbed by the termination. 


Conditions (1) and (2) can be illustrated by the extreme cases of termination of a line by an open 
circuit and a short circuit. Figure 7.5 shows what happens when there is an open circuit. Figure 7.5(a) 
shows the situation just before the incident wave has reached the open circuit. When it reaches it a 
reflected wave is generated to satisfy the boundary conditions. From Equation (7.37) the amplitude of 
the reflected wave is equal to that of the incident wave. A short time later the situation is as shown in 
Fig. 7.5(b), with the reflected wave superimposed on the incident wave. 
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Open circuit 


{b) 


Fig. 7.5 The reflection of a step wave by an open circuit. The wave is shown (a) just before it reaches 
the end of the line, and (b) just after it has been reflected. 
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In the same way Fig. 7.6(a) shows a wave approaching a short circuit. This time the wave is inverted 
on reflection so that the sum of the wave amplitudes is always zero at the short circuit. A short time 
after the reflection of the wave the situation is as shown in Fig. 7.6(b). 


Short circuit 
(b) 


Fig. 7.6 The reflection of a step wave by a short circuit. The wave is shown (a) just before it reaches 
the end of the line, and (b) just after it has been reflected. 


It is instructive to consider the effects of these reflections at the input of the transmission line. To 
make things simpler we will assume that the source is matched to the line so that Zs; = Zp as shown in 
Fig. 7.7(a). Figure 7.7(b) shows how the voltage at the start of the line varies with time when the 


termination is an open circuit. InitiallyV,, =V, /2as the source feeds current into the line to set up the 


incident wave. This wave travels down the line in time 7. At time 27 the reflected wave returns to A 


so thatV,, =V; . The reflected wave is completely absorbed by the source impedance, so no further 


change takes place and the line is charged to the full source voltage. 


P<) A 
V; R 
(a) 

V, Va 

V, v 
vs 2 

2 

t 
ft) t : 2T 
(b) aT {c) 


Fig. 7.7 Variation of the voltage at the input of a transmission line with time. The line has 
characteristic impedance Z,, and is connected to a matched source and to a mis-matched load. When 
the termination is an open circuit the voltage varies as shown in (b). (c) shows what happens when R 
is a short circuit. 
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If the line is terminated by a short circuit the voltage varies with time, as shown in Fig. 7.7(c). There 
is an initial period when the voltage at A is not zero, but this is cancelled at time 27 by the return of 
the reflected wave. According to elementary circuit theory the voltages would take their final values 
as soon as the switch was closed. The true situation can be described by saying that, until the return of 
the reflected wave, the source has no information about the magnitude of the terminating impedance. 
The step cannot travel faster than the speed of light, so it always takes a finite time for the system to 
reach a steady state. 


The transmission of a short pulse along a line can be investigated by considering the superposition of 
positive and negative step functions separated by a short interval of time. 


7.8 Transformation of impedance along a transmission line 

Consider a transmission line having characteristic impedance Zp which is terminated at x = 0 by an 
impedance Z, (see Fig. 7.8). At the point on the line which is / from the load, x = - / and the voltage 
is given by 


V =V,exp[j(or + k1)|+TV, exp| j(ot —k/)| (7.38) 


The corresponding expression for the current is 


r-( yi JeseLorsan]-r(Z JeseL (ont) (7.39) 

Zo Z, 

where 

pSee a (7.40) 
A ae 

from Equation (7.25). 


The apparent impedance at this point on the line is given by 


V _ exp( jk/)+T exp(—jkl) 
I exp(jkl)-T exp(—jdl) 


Z= Zi (7.41) 
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After substituting for l and rearranging, we get 


Zi _ Z,[exp(jkl)+ exp(- jk)]+ Z,[exp(jkl)- exp(- j)] 


ZZ, [exp( jk) —exp(— jkd)]+ Zo[exp( ski) + exp(- jk1)] 

_ 4 cos(kl)+ jZ, sin(k/) (7.42) 
Z, cos(kl)+ jZ, sin(k/) , 
Z, + jZ, tan(kd) 


Z, + jZ, tan(kd) 


This is a very important result. It shows that if the length / of the transmission line is of the same order 
of magnitude as the wavelength (2z/k), or greater, then the apparent impedance at the input of the line 
is not normally equal to the load impedance. The usual assumption of circuit theory is that the 
impedances are independent of the lengths of the connecting wires. We now see that this is just the 
limit of Equation (7.42) when A/ tends to zero. 


The significance of Equation (7.42) can be explained in a slightly different way by considering Fig. 
7.8. This shows two transmission lines, one of which is longer than the other by /. The longer line is 


terminated by a load Z, and the shorter by Z; , where Z;, is given by Equation (7.42). At the plane A—A 
the impedance presented to the incident wave is the same on both lines. Thus the reflected wave to the 
left of A-A must have the same amplitude and phase on both lines. It follows that Z, and Z; must 


reflect the incident wave with the same amplitude but with a phase difference of exp (2jkl), since that 
is the phase change from A-—A to Z, and back. 
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A 


Fig. 7.8 The apparent terminating impedance of a transmission line varies with position along the line. 
If two lines are terminated by Z,; and Z', as calculated from Equation (7.42), then the amplitude and 
phase of the reflected wave to the left of the plane A—A is the same on both lines. 


7.9 The quarter-wave transformer 


A case of considerable practical importance is obtained when the line is exactly one quarter of a 


wavelength long. Setting k/ = 2/2 in Equation (7.42) gives 


a 


pane Se 

Zo Z, 

or 

L227, (7.43) 


This result allows us to match a load to a line by inserting a short section of a line, having a different 
impedance, between them, as shown in Fig. 7.11. The short length of line must be a quarter 
wavelength long, and it must have the characteristic impedance given by (7.43). The same technique 
can be used to ensure a match between two transmission lines which have different impedances. 


| A/4 | 


Fig. 7.9 A source of impedance Z', can be matched to a load of impedance Z, by using a quarter 
wavelength of a transmission line whose characteristic impedance is Z) = V(Z, Z’,) as a transformer. 
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The quarter-wave transformer suffers from the difficulty that it works exactly only at a single 
frequency. This problem can be overcome by cascading a number of quarter-wave transformers to 
produce a broad-band match. 


Two special cases arise when the terminating impedance of the transformer is either an open circuit or 
a short circuit. If there is an open circuit, then Z, = 0 and (7.42) becomes 


Zw 4 


Zi, 7 jtan(k/) 


so that Z', tends to zero as k/ tends to 7/2. Thus a quarter-wave transformer transforms an open circuit 
into a short circuit. The converse is also true. If the termination of the transformer is a short circuit, 
then the impedance at its input tends to infinity. In either case there is a total reflection of the incident 
power. 


7.10 The coaxial line 


The simplest practical transmission line from the field point of view is the coaxial cable. We shall 
therefore take it as the example to be discussed, but the methods applied, and the conclusions reached, 
are applicable to all types of uniform two-wire line. A typical cable has a solid copper inner conductor 
and a braided outer conductor separated by a uniform cylindrical layer of polythene. To make the 
problem tractable we assume that the braided outer conductor can be represented with sufficient 
accuracy by a continuous conducting cylinder. Taking the radii of the inner and outer conductors as a 
and b, respectively, we recall that the capacitance per unit length is given by (6.50) 


Cc =— (7.44) 
and the inductance per unit length is given by (6.37) 


2a a 


L= fo ®) (7.45) 


Substituting these expressions for C and L into Equations (7.11) and (7.17) we get 


i (7.46) 


and 
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Z)= ~ J (2) (7.47) 
+h é a 


A very interesting result is obtained for the special case of an air-spaced line when 


>» = V Jeu 


1 
(8.854x10"? x 47x10)? 
0.2998 x10? ms 


= 
| 


II 


which is the experimental value for the velocity of light in free space. This suggests the existence of a 
link between electromagnetism and optics, and the possibility that light is an electromagnetic 
phenomenon. Since the relative permittivity of any medium is greater than unity, it follows that the 
phase velocity of the waves on a transmission line is always less than or equal to the velocity of light. 


7.11 The electric and magnetic fields in a coaxial line 


To find expressions for the electric and magnetic fields within a coaxial line which is carrying a wave 
along it we apply Gauss’ theorem and the magnetic circuit law. This approach involves us in making some 
assumptions about the fields in the line. Taking the electric field first, we assume that, as in the electrostatic 
case, only the radial component is present. The charge per unit length on the centre conductor is given by 


q=CV, exp| j(wt—/z)| (7.48) 
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so, applying Gauss’ theorem 


EB. = eae 1 expli(or—kz)] 
2mér 2nEr 


Substituting for the capacitance per unit length from Equation (7.44) gives 


eee ae 
7 aed ad kz)] (7.49) 


The magnetic field is given by 


I 1 V, 
H,=——=— + (at —k. 
oan Dae TZ exp| j(o—kz)| 


Substituting for Z) from Equation (7.47) 


E V; 1 é 
H,= /——* (@t—kz)\= |—E, 7.50 
=n ilejayr POEM VF = 


These fields are shown in Fig. 7.10. We note that the ratio of the strengths of the electric and magnetic 
fields is 


E 
Fi (ng (7.51) 
ily, & 


This quantity, which depends only upon the properties of the material filling the line, has the 
dimensions of resistance. It is known as the wave impedance of the wave. 


—— Electric field (£) 
“———™ Magnetic field (H,) 
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Fig. 7.10 The electromagnetic field in one wavelength of a coaxial line showing the directions of the 
electric and magnetic fields and the direction of propagation of the wave. Half of the outer conductor 
has been removed for clarity. The electric field lines are shown schematically; they actually radiate 
equally in all directions. 


Although these results have been derived by making assumptions about the directions of the fields, 
they are confirmed by more rigorous analysis. It is found that, for any two-wire transmission line, the 


electric and magnetic field vectors are perpendicular to each other and to the direction of propagation 
of the wave. Such waves are known as transverse electric and magnetic (TEM) waves. 


7.12 Power flow in a coaxial line 
When discussing the storage of energy in capacitors and inductors we saw that it is sometimes useful 
to think of the energy as being distributed throughout the electric or magnetic fields. In the same way 


we can think of the power flowing down a coaxial cable as being distributed throughout the field. At 
any point the electric energy density is 


1 
We =ZD-E=7ek, (7.52) 


while the magnetic energy density is 


We = 5 BH = xt H; (7.53) 


Making use of Equation (7.51) we can write the total energy density as 


1 
W=Ww,twWy = 5 eB +om [£e7|=ee: (7.54) 
0 


Now at any point within the line the fields vary sinusoidally with time. Therefore the time average of 
the stored energy is 


a ee EVe 1 
eS |, es 7.55 
w= 968 = Tina) er 


The whole field pattern is moving in the positive z direction with velocity v,, so the power density is 


1 1 
S=Wv, =5E, (ev, E,)= se Hy (7.56) 


137 


Electromagnetism for Electronic Engineers 7. Transmission lines 


making use of Equations (7.46) and (7.51). We shall see in Chapter 8 that this is an example of the 
application of a more general result relating to the power flow in electromagnetic waves. 


To check the equivalence of the field and circuit descriptions of the problem we can calculate the total 
power flow by integrating Equation (7.56) over the space between the conductors. Then 


dr (7.57) 


{’ ev Vo Arr 
@2(In(b/a)) rr? 


2 
p= ye {= 2x Ns (7.58) 
2 Ly In(b/a) 2Z, 


from the expression for Zp given in Equation (7.47), so the power is identical to that calculated by 
circuit methods (Equation (7.27)). 


7.13 Summary 


In this chapter transmission lines have been introduced through a discussion of distributed circuits. 
The voltages and currents were found to be governed by the wave equation, and it was demonstrated 
that solutions could be found which could be interpreted as travelling waves. The power carried by 
these waves was expressed in terms of the voltage on the line and its characteristic impedance. 
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It was demonstrated that, unless a line is matched by being terminated by a load equal to its 
characteristic impedance, some of the power in a wave incident upon the termination is reflected back 
down the line. This reflection produces a partial standing wave on the line. The effect of the reflected 
wave is to make the input impedance of a line depend upon its length as well as upon the terminating 
impedance. Thus it was shown that, at high frequencies when the lengths of lines are comparable with 
the wavelengths of the signals on them, it is no longer possible to treat circuits in terms of lumped 
components connected together by wires of arbitrary length. The transformation of impedance by a 
short length of a transmission line is put to use in the quarter wave transformer to match a load to a 
line. 


It is also important for lines to be correctly matched when the signals are in the form of short pulses. 
Incorrect matching results in reflections of the pulses. If a line is not matched at both ends multiple 
reflections can occur so that a whole train of pulses is produced by transmitting a single pulse down 
the line. 


Transmission lines were also considered from the point of view of the electric and magnetic fields 
within them. It was shown that the phase velocity and characteristic impedance of a line at a given 
frequency can be calculated from field considerations. Expressions were found for the electric and 
magnetic fields within a coaxial line, and it was demonstrated that the flow of power down the line 
could be regarded as that of an electromagnetic wave propagating in the space between the 
conductors. It was found that, for an air-spaced line, this wave propagates with a phase velocity equal 
to the speed of light. The equivalence between the circuit and field approaches to the description of 
waves on transmission lines was demonstrated by showing that both give the same expression for the 
power flow. 
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8. Maxwell’s equations and electromagnetic 
waves 


Objectives 

e To show how Maxwell removed an inconsistency in the magnetic circuit law by introducing 
the idea of the displacement current. 

e To derive the differential forms of the magnetic circuit law and Faraday’s law, and to 
introduce the curl of a vector. 

e To present Maxwell’s equations in both integral and differential form. 

e To demonstrate how the existence of plane electromagnetic waves in free space can be 
deduced from Maxwell’s equations. 

e To show that the power flow in an electromagnetic field can be calculated by integrating the 
Poynting vector over a closed surface. 


8.1 Introduction 


In Chapter 7 we saw how it is possible for electromagnetic waves to propagate along a transmission 
line such as a coaxial cable. We shall now proceed to show that electromagnetic waves can exist 
independently of any system of conductors and can travel through free space. This conclusion was 
first published by James Clerk Maxwell in 1873. The set of equations which bears his name forms a 
summary of all the topics treated in this book and demonstrates the symmetry and unity of 
electromagnetic theory. At the same time they are the foundation upon which modern electromagnetic 
theory is built. We are now so accustomed to the presence of radio, television, mobile phones and 
radar in the world that it is hard to realize how great Maxwell’s contribution was to the understanding 
of the theory of electromagnetism. In his own day his theory was regarded as quite outrageous by 
many until its validity was demonstrated by Hertz’s experiments some fifteen years later. 


8.2 Maxwell’s form of the magnetic circuit law 


Maxwell noticed that, under certain circumstances, the magnetic circuit 
law in the form derived in Chapter 4 can give inconsistent results. This 
is readily demonstrated by considering the circuit shown in the figure on 
the right — a charged capacitor discharging though a resistor. 


To apply the magnetic circuit law we must define a closed path C 
which encircles the wire and an open surface S which spans that Surface S 
path as a soap film can span a loop of wire. The figure in the 
margin attempts to show what is meant. If S'is chosen so that it 
passes through the wire, then no problem arises and the magnetic 


field satisfies Equation (5.7), which is repeated here for Path C 


convenience: 
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f H-di={{J-da (8.1) 
Cc S 


But, if the surface S is chosen so that it passes between the plates of the capacitor, then J = 0, 
implying that the line integral of H around C is also zero. This sounds nonsense! We expect the 
magnetic field around the circuit to be independent of the details of the means by which it has been 
calculated. 


To resolve this difficulty we consider two different surfaces S;and S; which both span C, as shown in 
Fig. 8.1, and which together form a closed surface. 


I, ¥ 


S; 


t, ¥, 


Fig. 8.1 Maxwell showed that the conduction current (/) passing through a closed loop (C) can be 
dependent upon the choice of the surface (S) over which it is calculated. The sum of the conduction 
and displacement currents (/ + d¥/dt) is independent of the choice of surface. 
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At a particular instant let the total charge enclosed within this surface be Q, and let currents I, and /, 
cross the two surfaces as shown. Now charge must be conserved, so the rate of change with time of 
the charge enclosed must be equal to the net current flow into the volume enclosed. That is, 


d 
Cais “ (8.2) 


If we define ¥Y; to be the flux of D through S, and ¥ that through S,as shown in Fig. 8.1, then the 
application of Gauss' theorem gives 


Y -Y,=2 (8.3) 


(8.4) 


This equation shows that the quantity (7 +da¥/ dt) is independent of the choice of the position of the 


surface S even when / is not. This quantity is known as the total current, and d‘¥/dt is known as the 


displacement current. 


When the current density on the right-hand side of Equation (8.1) is replaced by the total current 
density (J + dD/dt), the result is the equation 


dD 
frt-al=[j{s+)-aa (8.5) 


This form of the magnetic circuit law is universally true. 


8.3 The differential form of the magnetic circuit law 


In Chapter | we saw that it is possible to express Gauss' theorem in either an integral form (Equation 
(1.5)) or a differential form (Equation (1.9)). The differential form turns out to be the more useful of 
the two because it enables a wider range of problems to be solved. The same is true of the magnetic 
circuit law, as we shall see. 


To derive the differential form of Equation (8.5) we apply it to the small rectangular path ABCD 
shown in Fig. 8.2. The surface S is taken to be in the y—z plane. The line integral of H is taken around 
ABCD in a direction which is in the right-handed corkscrew sense with respect to the positive x- 
direction. In other words, the elementary loop is encircled in the direction ABCD. 
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x 


Fig. 8.2 The differential form of the magnetic circuit law can be derived by considering the line integral 
of the magnetic field around the infinitesimal loop ABCD. 
The value of the line integral of H around the loop is computed by considering each side in turn. The 


effect of the dot product of H with dl is to pick out the component of H which is parallel to the path. 
Thus, for AB we get 


[ Heal =H, 5 (8.6) 


Along CD the y-component of H can be written (7 at (oH } dz)&), so that 


[. Hat-| {Zlalem (8.7) 


The minus sign on the right-hand side of this equation is there because the path CD is traversed in the 
negative y-direction. 


Combining Equations (8.6) and (8.7) gives 


B D 1- OH, 
| Hdl = 1. H-dl=- ; VEZ (8.8) 


IZ 


Applying the same argument to the other two sides of the loop gives 
Cc A OH 

| H-dl + \ H-dl=—. 5v& (8.9) 
B D oy 


The complete line integral around the loop ABCD is found by adding Equations (8.8) and (8.9) to give 


HH 
f nat=| Ze Fla (8.10) 
Oy Oz 
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Since the loop is small, it is reasonable to assume that J and D are constant over it, so that the integral 
on the right-hand side of Equation (8.5) becomes 


aD aD, 
— |-dA=| J, +—= 8.11 
IIs >). [v4 a) ava (8.11) 


This time the dot product picks out the components of J and D which are normal to the plane of the 
loop because the vector area dA is also normal to the plane of the loop. 


Equating Equations (8.10) and (8.11) yields the differential form of Equation (8.5) for the x-direction: 


_ eo (8.12) 


Similar equations can be obtained for the other two coordinate directions: 


oD 
a Aas (8.13) 
Oz Ox ' Ot 
OH 
y OH, _, , OD. (8.14) 
Ox Oy "Ot 
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Equations (8.12) to (8.14) together form the differential version of Equation (8.5) expressed in terms 
of the components of the vectors. 


As usual, it is possible to write these equations in a more compact form by using vector notation. To 
do this we multiply each by the appropriate unit vector and add them together to give 


OH, OH,), (OH, OH,\. (OH, OH,), oD 
an x+ 2S | y+ —-—= |z=J+—— (8.15) 
Oy Oz Oz Ox Ox Oy 


The left-hand side of this equation is known as the curl of H and can be written as a determinant: 


ky a 
curl H = Bs aca ca (8.16) 
Ox Oy 
iW. i. HH, 
which can be recognized as the vector product of the vector operator V and the vector H. Thus 
Equation (8.15) can be written succinctly in the form 
oD 
VAH=J+— (8.17) 
Ot 
Like the other differential forms derived earlier, Equation (8.17) is actually valid for coordinate 
systems other than the rectangular Cartesian system used to derive It. 
8.4 The differential form of Faraday’s law 
The general integral form of Faraday’s law was shown to be 
OB 
E-dl=-|| —-dA (8.18) 
es es 


(Equation (6.18)) for the case when the circuit is fixed in space and the magnetic flux density is 
changing with time. 


So far we have assumed that the line integral of E is to be taken around a loop of wire. But if the loop is 
open-circuited and the potential difference between its ends is measured with a. high-impedance volt- 
meter, then the loop can be thought of as a device for measuring the value of the line integral of an electric 
field which exists in space whether the wire is present or not. Thus Equation (8.18) can be regarded as a 
generalization of Faraday’s law, implying that, if there is a changing magnetic field in any region of space, 
then there is also an electric field there. This is, strictly speaking, a plausible guess rather than a deduction, 
but it is one whose validity has been demonstrated by the correctness of the results derived from it. 
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Making a comparison between Equations (8.5) and (8.18) allows us to deduce straight away that the 
differential form of (8.18) must be 


Vans o” (8.19) 
Ot 


8.5 Maxwell’s equations 


We have now introduced all the principles of electromagnetism and it is convenient to gather the main 
results together. Taking the integral forms first they are: 


Gauss' theorem in electrostatics (Equation (2.5)) 

ff D-ds=|[[ pav (8.20) 
Gauss’ theorem in magnetostatics (Equation (4.15)) 

{} B-dA =0 (8.21) 


The magnetic circuit law (Equation (8.5)) 


dD 
fueal=ff{s+2)aa (8.22) 


Faraday’s law of induction (Equation (6.18)) 


eB 
pE-di=-[[—-da (8.23) 


The corresponding differential forms are 


(2.6) V-D=p (8.24) 
(4.16) V-B=0 (8.25) 
(8.17) VaH=J+< (8.26) 
(8.19) VAE --= (8.27) 
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These four equations form a summary of the whole of fundamental electro-magnetic theory. They are 
known, collectively, as Maxwell’s equations and are the starting point for the discussion of all the 
more advanced topics in electro-magnetism. 


In order to make use of them we also need a number of other equations which have been introduced in 
previous chapters. First there is the continuity equation 


(3.14) yy Je (8.28) 


Then there are the constitutive relations which introduce the properties of materials, albeit in an 
idealized form. These are: 


(2.4) D=cE (8.29) 
(3.3) J=oE (8.30) 
(5.5) B=uH (8.31) 


The application of Equations (8.24) to (8.31) is most easily accomplished by making use of the 
methods of vector calculus. These techniques lie beyond the scope of this book, but a simple example 
of the use of Maxwell’s equations which can be dealt with using more elementary mathematical 
methods is the subject of the next section. 
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8.6 Plane electromagnetic waves in free space 


In the previous chapter we saw that electromagnetic waves can propagate along a coaxial line and that 
the electric field vector, the magnetic field vector, and the direction of propagation are all mutually 
perpendicular. Taking our clue from this, we can see whether Maxwell’s equations indicate that 
similar waves can propagate in free space far from any material boundaries. 


Let us assume that the wave propagates in the z-direction, and that the electric field has a component 
only in the x-direction. We assume, furthermore, that the intensity of the electric field varies only in 

the z-direction. Then from Equation (8.27), making use of the definition of the curl of a vector given 
in Equation (8.16), we have 


4 == (8.32) 


since all the other components of V x E are zero. This equation shows that the magnetic field vector 
must lie in the y-direction, at right angles to both the electric field and the direction of propagation. 
Thus Equation (8.32) can be written as a scalar equation 


OE. OB, 
ee (8.33) 
Oz Ot 
In free space the conduction current J = 0, so Equation (8.26) becomes 
OH, OD, 
—-— =—__ (8.34) 
Oz Ot 


The magnetic field terms can be eliminated between these two equations by differentiating the first 
with respect to z and the second with respect to ¢, and by making use of Equations (8.29) and (8.31) 


with €=€&, and “= 4). The result is that 


= = ee) (8.35) 
<= =6§ ; ; 
aS Beet BP 
That is 
OE 10°E 
x = x 8.36 
Oz? Vv. Ot? ( ) 
where 
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1 
yo= (8.37) 


° Vv e0Ho 


is the velocity of light. Equation (8.36) is the wave equation which was first encountered in Equation 
(7.4). 


If the electric field terms are eliminated from Equations (8.33) and (8.34) in a similar manner, the 
result is 


ee (8.38) 


The two wave equations have the following general solutions for waves travelling in the positive z- 


direction: 
E. = E, exp| j(ot —kz)| (8.39) 
and 

H, =H, exp j(wt-kz)| (8.40) 


where Ep and Hp are complex amplitudes which incorporate the relative phases of the electric and 
magnetic fields. 


The relationship between the two fields can be found by substituting Equations (8.39) and (8.40) into 
either Equation (8.33) or Equation (8.44). Making use of Equation (8.33), we get 


— jkE,, exp| j(at —kz)|=—jau,H, exp| j(ot —kz)| (8.41) 


From this it can be seen that £,and H, are in phase with one another. Also that 


Ey _ [Ho 3770 (8.42) 
A, €o 


We recognise this as the wave impedance of a wave in free space (see Equation (7.51)). It is 
sometimes referred to as the intrinsic impedance of free space. 


Figure 8.3 illustrates the complete solution which has been obtained. It is necessary to be a little 
careful in interpreting this diagram. It shows the amplitudes of the field vectors rather than being a 
map of the field. The whole field pattern is moving in the z-direction with the velocity of light. At any 
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instant the electric and magnetic field strengths are uniform at all points in a plane perpendicular to 
the z-axis. 


Fig. 8.3 The electric and magnetic fields in a plane electromagnetic wave propagating in the z- 
direction. 


150 


Electromagnetism for Electronic Engineers 8. Maxwell’s equations and electromagnetic waves 


8.7 Power flow in an electromagnetic wave 


In Chapter 7 it was shown that the power flow in a coaxial line could be calculated correctly by 
assuming that the stored energy in the electromagnetic field travelled with the phase velocity in the 
direction of the wave. For waves in free space the energy density is, from Equation (7.54), 


w=e,E? (8.43) 


The instantaneous power density in the wave is 


=F H, (8.44) 


Now S is in the z-direction, so it is possible to write Equation (8.44) in vector notation as 
S=EAH (8.45) 
The vector S is known as the Poynting vector. Poynting’s theorem states: 


The integral of (E A H) over a closed surface is equal to the instantaneous flow of electromagnetic 


power out of the volume enclosed by that surface. 


The proof of this theorem is beyond the scope of this book. Note that, strictly speaking, it is the 
integral of S over a closed surface which has meaning rather than S itself. We have already seen that 
the integral of the Poynting vector gives the correct answer for the power flow in a lossless coaxial 
line. 


8.9 Summary 


We have considered how Maxwell removed an inconsistency in the magnetic circuit law by 
introducing the idea of the displacement current. The modified law was expressed in differential form 
by introducing the curl of the vector H. Faraday’s law was also put into differential form by using the 
curl of E. 


The mathematical statements of the laws of electromagnetism were collected together in both their 
integral and differential forms. This set of equations is known as Maxwell’s equations. It was 
demonstrated that, for the special case of plane waves in free space, Maxwell’s equations lead 
naturally to the prediction of the existence of electromagnetic waves which travel with the speed of 
light. 


Finally, it was shown that the power flow in the electromagnetic field can be represented by the 


Poynting vector (E A H) : 
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Appendix 


Physical constants 

Primary electric constant (9) 
Primary magnetic constant (10) 
Velocity of light in vacuum 
Wave impedance of free space 
Charge on the electron 

Rest mass of the electron 
Charge/mass ratio of the electron 


Properties of dielectric materials 


$854: 10°" 
An x 107 
0.2998 x 10° 
376.7 

-1.602 x 10°? 
9.108 x 10°! 
-1.759 x 10"! 


Relative permittivity 


Alumina 99.5% 10 
Alumina 96% 9 
Barium titanate 1200 
Beryllia 6.6 
Epoxy resin 3.5 
Ferrites 13-16 
Fused quartz 3.8 
GaAs (high resistivity) 13 
Nylon 3.1 
Paraffin wax 2.25 
Perspex 2.6 
Polystrene 2.54 
Polystyrene foam 1.05 
Polythene 2.25 
PTFE (Teflon) 2.08 


Properties of conductors 


Conductivity (S m’') 


Aluminium 3.5 x 10’ 
Brass 1.1 x 107 
Copper 5.7 x 10' 
Distilled water lei 
Ferrite (typical) 10° 

Fresh water 10° 

Gold 4.1.x 10° 
Iron 0.97 x 107 
Nickel 1.28 x 10’ 
Sea water 4 

Silver 6.1 x 10’ 
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Steel 0.57 x 10’ 


Properties of ferromagnetic materials 


Ly Saturation magnetism (Bsa) (T) 
Feroxcube 3 1500 0.2 
Mild steel 2000 1.4 
Mumetal 80 000 0.8 
Nickel 600 
Silicon iron 7 000 1.3 


Summary of vector formulae in Cartesian coordinates 
a=a,X+ay+aL 
a-b=a,b, +a,b, +a,b, 
aAb=(a,b, —a_b, ie +(a,b, -—a,b, + (a,b, —a,b, hi 
Ox On. 0% 
V=—x+—yt+—Zz 
Ox Oy oz 
OV. OV OV 


rad V =VV =—x+—y+—z 
= Ox ay" Oz 
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diva=V-a=—*+ 


da, Oa,|, | ey Oa, @a,|, 
curla = V Aa=4—*—-—— -K +5 —* -—* +y +) -—* 71 
Oy oz Oz Ox Ox oy 


OV AV av 
= + + 
ax? Gy? Gz” 


VV 


Summary of the principal formulae of electromagnetism 
Inverse square law of electrostatic force 


F = 22, r 
ATE r 
Relationship between E and D 
D=cE=6,¢,E 
Gauss’ theorem 


ff D-ds= |] pdv divD=V-D=p 


Electrostatic potential difference 
B 
V,-V,=—-( E-dl E =~ gradV =-VV 


Poisson’s equation 
i 
é 
Energy density in an electric field 
w= a D-E 
2 
Ohm’s law 
J=oE E=pJ 
Power dissipated per unit volume 
p=E-J 


Continuity equation 
ff s-aa=-£ fff pav Vaiee! 
S Oty 


Electromotive force 
~=fEdl 


Law of force between moving charges 


Fe 20, Pa MyQQ (v, a(v, a8) 
Are r™ Arr 
Force on a moving charge 


F = O(E+(vB)) 
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Force on a current-carrying conductor 
F = (dl ,B) 


Biot-Savart law 


Mol ¢ dlar 
B=— 
Ar i r 
Relationship between B and H 
B= WH = 1 ,H 
Magnetic circuit law as modified by Maxwell 
D 
fueal=ff{s+)-aa cwlH=VaH=5+22 
C Ss dt Ot 
Conservation of magnetic flux 
{} B-dA=0 divB=V-B=0 
Faraday’s law of induction 
B 
fe-d=—|fB-aa cwlE=VAE=-2 
Ot Ot 


Energy density in a magnetic field 
| 

w=—B-H 
2 

The Poynting vector 

S=EAH 
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